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PREFACE 


This  work  was  conducted  under  sponsorship  of  the  Federal  Aviation 
Administration  (FAA)  as  part  of  the  runway  roughness  study  to  determine 
the  effect  of  the  loaded  pavement  profile  on  aircraft  operations.  An 
outgrowth  of  this  study  is  a noncontact  system  to  measure  the  pavement 
deflection  basin  under  a moving  load.  This  report  describes  the  hard- 
ware developed  and  associated  transfer  function  theory  which  forms  the 
basis  for  a pavement  evaluation  and  response  scheme.  (The  effect  of  the 
loaded  pavement  profile  will  be  the  subject  of  a separate  report  to  be 
published  in  September  1977.)  Continued  development  of  the  noncontact 
system  reported  herein  is  being  sponsored  by  the  FAA  under  separate 
contract. 
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Today,  personnel  concerned  with  the  evaluation  of  airfield 
pavements  are  confronted  with  three  major  problems.  First,  many 
primary  pavements  are  old  and/or  are  approaching  early  stages  of 
deterioration.  Second,  the  gross  weight  of  aircraft  is  increasing 
continually,  thus  the  demands  on  airfield  pavements  are  also 
increasing.  Third,  procedures  (References  1,  2,  3 and  4)  used  to 
evaluate  airfield  pavements  are;  either,  destructive  in  nature,  or, 
apply  to  only  small  areas  of  the  pavement:  in  addition,  they  require 
considerable  performance  time  and  are  therefore  costly  with  respect 
to  the  flow  of  traffic. 

The  problem  faced  by  the  Air  Force  in  this  context  has  been 
stated  by  one  author  (Reference  5)  who  wrote, 

"The  Air  Force  alone  owns  enough  pavement  to  be  able 
to  provide  a 200-foot  wide  runway  stretching  from  the 
state  of  Washington  to  the  Southern  tip  of  Florida.  The 
problem  of  maintenance  and  rehabilitation  of  this 
inventory  becomes  more  complex  and  critical  each  year. 

Most  of  the  systems  are  over  20  years  old  . . . coupled 
with  these  aging  pavements  is  the  rapid  growth  of 
aircraft  traffic  and  weights." 

The  need  exists  today  for  a method  of  pavement  evaluation  which 
would  be  able  to  quantify  the  stage  of  the  aging  process  of  a pavement 
and  provide  guidance  with  respect  to  its  rehabilitation.  It  should  be 
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able  to  do  so  considering  the  action  of  varying  magnitudes  and 
configurations  of  loads,  and  under  varying  ambient  conditions. 

Ideally,  it  should  be  nondestructive,  rapid,  permit  the  evaluation 
of  the  entire  pavement  with  a minimum  interruption  to  air  traffic, 
and  should  be  inexpensive  and  easy  to  use. 

The  problem  of  evaluating  a pavement  is  a complex  one.  The 
pavement  section  consists  of  various  materials.  These  materials  are 
far  from  the  ideal  models  of  classical  mechanics ; and  they  vary 
diumally,  seasonally  and  with  repetitions  of  loading.  In  addition, 
both  vehicular  and  non~vehicular  loads  applied  to  the  pavement- 
subgrade  system  vary  in  magnitude,  intensity,  and  frequency. 

Recognizing  the  nature  and  scope  of  the  problem,  and  the  need 
that  exists  today  to  overcome  it,  research  activities  were  initiated 
at  Purdue  University  toward  the  development  of  a nondestructive 
pavement  evaluation  capability  to  account  for  the  complexity  and 
variability  of  the  pavement-subgrade  system,  and  simultaneously  satisfy 
the  demands  imposed  by  practical  considerations;  such  as  mobility  and 
ease  of  operation  of  equipment.  Research  in  this  endeavor  is  predicated 
on  the  use  of  transfer  function  theory  as  the  key  to  pavement 
evaluation.  This  methodology  has  evolved  over  a 15  year  period  from 
concept,  through  theory  and  laboratory  studies  to  field  investigations 
using  fixed  installations  (References  6,  7,  8,  9,  10,  11,  12  and  13). 


Yards  and  Docks  (Reference  2);  the  Canadian  Department  of  Transporta- 
tion (Reference  3);  the  British  Department  of  the  Environment 
(Reference  4)  and  the  U.  S.  Department  of  the  Army  and  the  Air  Force 
(Reference  18).  Basically,  the  five  procedures  involve  the  use  of 
field  tests  in  conjunction  with  sampling  and  laboratory  testing. 

The  FAA  procedure  uses  California  Bearing  Ratio  (CBR)  tests 
conducted  in  accordance  with  procedures  in  MIL- STD  621A  Method  101. 
The  Navy  method  uses  non-repetitlve  plate  load  tests  conducted  accord- 
ing to  specifications  of  American  Society  for  Testing  and  Materials 
(ASTM)  Test  Method-Designation  D-1196  (Reference  19).  The  Canadian 
practice  uses  repetitive  plate  bearing  tests  specified  in  ASTM 
Test  Method  Designation  D-1195  (Reference  20).  The  British  proced- 
ure also  uses  repetitive  plate  bearing  tests  but  these  are  conducted 
in  a slightly  different  manner  from  the  ASTM  D-1195  Test  Method. 

The  Air  Force  evaluation  procedure  is  summarized  in  the  report  of 
the  testing  program  conducted  in  a recent  (1975)  evaluation.  It  was 
reported  (Reference  21)  that, 

"Field  testing  consisted  of  28  exploratory  test 
pits  . . „ Thickness  measurements  were  made  on  each 
pavement  component,  in  situ  density  and  moisture 
contents  were  determined,  values  of  Modulus  of  Sub- 
grade Reaction  and  CBRs  were  obtained  and  bulk 
samples  of  soil  layers  were  taken.  Core  samples 
were  extracted  at  184  locations  to  ascertain 
pavement  thickness  and  subsurface  profile.  Lab- 
oratory testing  included  classification  of  soils; 
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development  of  soil  moisture-density  relationships; 

Marshall  testing  of  asphaltic  concrete;  and  tensile 
splitting  of  Portland  cement  concrete  cores." 

The  period  of  the  evaluation  was  about  8 days. 

The  prevalent  feeling  concerning  the  yse  of  test  pits,  and 
current  destructive  evaluation  procedures  was  quite  aptly  put  when 
one  author  (Reference  22)  wrote, 

"Yet  most  civil  engineers  still  cannot  scientifi- 
cally evaluate  the  load  carrying  capacity  of  an 
airport  pavement  without  first  rendering  it  unservice- 
able in  a way  no  aircraft  could  have  accomplished. 

The  resulting  test  pits  when  filled  and  resurfaced 
remain  as  abrupt  discontinuities  to  uniform  pavement 
performance." 

This  general  feeling  coupled  with  the  many  other  disadvantages  of 
destructive  testing  has  directed  research  towards  the  development 
of  rapid,  nondestructive  evaluation  equipment  and  procedures. 

2.1.2  Nondestructive  Evaluation 

Nondestructive  procedures  are  based  on  the  analysis  of  the 
measured  surface  deflection  response  of  a pavement-subgrade  system 
to  loads  applied  at  the  surface.  The  main  differences  among  existing 
methods  are:  (1)  the  magnitude  and  manner  in  which  the  load  is 
applied  to  the  surface,  and  (2)  the  particular  aspect  of  the  surface 
disturbance  that  is  measured,  i.e.,  deflection  magnitude,  wave  length, 
frequency.  Table  1 shows  a summary  of  the  major  load  classes  and  their 
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TABLE  1 SUMMARY  OF  NONDESTRUCTIVE  TESTING  EQUIPMENT  AND 
PERTINENT  CHARACTERISTICS 


Category 

Load  Class 

Device 

Deflection 
Measurement  Made 

1 

Static 

Plate  (19)a 

Magnitude 

Plate  (20) 

Magnitude 

2 

Impact 

Washington  State  (58) 

Magnitude 

Paris  (58) 

Magnitude 

German  (58) 

Magnitude 

3 

Vibratory 

Shell  (23) 

Wave  Length 

Road  Rater  (25) 

Magnitude 

Dynaflect  (23) 

Magnitude 

Corp  of  Engr.  (25) 

Magnitude 

AFCEL  (68) 

Wave  Length 

4 

Vehicular 

Benhleman  (69) 

Magnitude 

California  (55) 

Magnitude 

Lacroix  (70) 

Magnitude 

TRRL  (71) 

Magnitude 

South  Africa  (23) 

Curvature 

lumbers  refer  to  references 
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characteristics,  devices  being  developed  and  In  use,  and  the 
particular  aspect  of  the  surface  disturbance  that  is  measured. 

In  the  static  load  class  repetitive  or  non-repetitive  loads 
are  applied  to  plates  of  varying  sizes  on  the  surface  of  the  pavement. 

In  the  Impact  load  class,  weights  are  allowed  to  free-fall  and  Impact 

a plate  on  the  pavement  surface.  The  weights  may  range  in  size  from 

■» 

20  lbs  (Reference  23)  to  500  lbs , as  was  used  by  Isaada  (Reference  24) . 

The  vibratory  procedures  also  show  a considerable  range  in  magnitude  of 
vibratory  loads  and  frequency  of  vibration.  Peak  vibratory  loads  may 
range  from  0.75  to  25  kips  and  frequencies  may  range  from  1 to  5000 
hertz.  Green  and  Hall  (Reference  25)  conducted  a thorough  study  of 
vibratory  equipment.  Their  work  and  the  references  cited  In  Table  1 
provide  additional  details.  Of  the  devices  that  use  vehicular  loads, 
the  Benkleman  beam  measures  rebound  deflection,  and  the  LaCrolx  and 
California  deflectometer  are  somewhat  automated  Benkleman  beams.  The 
TRRL  device  employs  an  optical  (non  contact)  displacement  transducer 
for  measuring  deflections,  and  the  South  African  apparatus  measures 

curvature  of  the  pavement  as  the  wheel  approaches.  j 

i 

2.2  Analysis  | 

Numerous  analytical  solutions  exist,  employing  various 

models,  representing  a pavement-subgrade  system.  These  models  are 

I 

based  on  the  theories  of  elasticity  and  viscoelasticity.  Numerous 

I 

equations  representing  their  behavior • under  static  and  dynamic  loading 
have  also  been  solved  (References  26  through  40).  In  1972,  Boyer  (Ref- 
erence 11)  presented  a rather  definitive  review  of  these  solutions.  In 
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addition,  he  employed  the  concept  of  transfer  function  theory  along 
with  the  results  of  full  scale  field  testing,  using  gages  installed  in 
the  pavement.  He  proposed  an  "a  posteriori"  modelling  concept  to 
develop  transfer  functions  and  predict  pavement  response.  Highter 
(Reference  13)  reviewed  the  use  of  energy  concepts  and  methods,  and 
failure  criteria  relative  to  pavement  evaluation.  Using  energy 
concepts  he  was  able  to  substantiate  that  "a  functional  relationship 
exists  between  cumulative  energy  as  measured  by  cumulative  peak 
deflections  imparted  to  a given  pavement  system  and  the  condition  of 
the  system".  Three  sources  of  data  were  studied  to  verify  his 
hypothesis.  AASHO  Road  Test  data,  the  Air  Force  base  traffic  records 
from  Pease  and  Castle  Air  Force  Bases,  and  data  from  field  testing  of 
overlays  specifically  constructed  for  his  investigation. 

A theoretical  study  of  dynamic  stiffness  and  its  application  to 
vibratory  nondestructive  testing  of  pavements  was  conducted  by  Weiss 
(Reference  41)  in  conjunction  with  the  study  by  Green  and  Hall 
(Reference  25) . Weiss  developed  a nonlinear  vibration  theory  for 
pavements  and  gave  a method  of  obtaining  the  shear  modal  us  and 
thickness  of  each  pavement  layer.  In  this  work  a vibrator  is  used  to 
generate  the  deflections  from  which  the  dynamic  stiffness  is 
calculated. 

Wiseman  (Reference  42)  used  the  Hogg  and  the  Hertz  models  to 
estimate  pavement  characteristics.  The  Hogg  model  provided  the 
required  pavement  stiffness  for  full-depth  asphalt  airfield  pavements. 
He  then  related  the  computed  stiffness  to  projected  traffic  and  a 
characteristic  length  using  published  Asphalt  Institute  procedures. 
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Wiseman  also  used  the  Hertz  model  to  compute  the  single  wheel  load  on  a 
rigid  pavement  as  a function  of  measured  stiffness. 

A probabilistic  analysis  of  the  response  of  pavement  structures 
was  undertaken  by  McCullough  (Reference  17).  He  cited  the  application 
of  stochastic  process,  in  this  regard  by  the  Texas  Highway  Department 
(Reference  43),  California  Division  of  Highways  (Reference  44)  and  the 
Asphalt  Institute  (Reference  45).  He  also  illustrated  the  application 
of  the  Markov  process  and  its  theoretical  relation  to  states  of  a 
pavement.  Reliability  techniques  have  also  been  reported  by  Darter 
and  Hudson  (Reference  46).  In  1976,  Moavenzadeh  (Reference  47) 
presented  a stochastic  simulation  model  for  predicting  pavement 
performance.  Recognizing  the  limitations  of  his  model,  he  concluded 
that  it  can  be  used  to  compare  various  design  alternatives.  It  also 
accounts  for  geometry,  load  area  and  Intensity,  random  inter-arrival 
times,  statistical  variation  in  material  properties  and  temperature 
histories. 

Inherent  in  all  analytical  and  stochastic  models  is  the  notion 
that  parameters  can  be  obtained  that  reflect  the  actual  vehicular  load 
experienced  by  a pavement  subgrade  system.  At  the  present  time,  the 
literature  showe  that  parameters  used  in  analysis  are  obtained  by 
static  field  test,  laboratory  testing,  vibrations,  impact  and 
measurement  of  pavement  rebound  characteristics.  The  actual  time 
dependent  response  of  the  pavement-subgrade  system  subject  to  one 
complete  vehicular  pass  is  not  considered. 
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SECTION  III 

THEORY  AND  METHODOLOGY 

There  exist  today  numerous  elastic  and  viscoelastic  solutions 
(References  26  through  40)  to  the  overall  vehicle-pavement-subgrade 
Interaction  problem,  in  addition  to  the  results  of  several  full  scale 
field  testing  programs  (References  48  through  50) . Accumulated 
theory  and  prototype  testing  point  to  the  following  fundamental 
observations : 

1.  The  characteristics  of  the  component  materials  of  a 
pavement-subgrade  system  vary  locally,  even  within 
the  small  volume  of  a representative  "homogeneous" 
sample . 

2.  Pavement-subgrade  characteristics  and  response  to 
loadings  change  with  seasonal  and  ambient  conditions. 

3.  Load  repetitions  affect  the  various  components  of  the 
system  differently. 

4.  A nondestructive  field  test  conducted  in  situ  using 
prototype  loading  better  reflects  the  action  of  a 
pavement-sub grade  system  than  destructive  tests  on 
samples,  or  tests  conducted  in  situ  using  non- 
representative  loads. 

5.  The  pavement  subgrade  system  exhibits  three  basic 
behavioral  characteristics;  namely, 

(a)  inertial  characteristics, 


(b)  elastic  characteristics,  i.e.  linear  and 


independent  of  time,  and 

(c)  viscoelastic  characteristics  that  account  for 
the  memory  and  the  time  dependent  nature  of 
the  system  response. 

Recognition  of  the  complexity  and  variability  of  the  pavement - 
subgrade  system  and  the  nature  of  the  wheel  loading-pavement-subgrade 
interaction  problem,  has  led  to  the  development  of  global  methods  of 
characterizing  the  pavement-subgrade  system  under  prototype  conditions 
of  loadings  (References  38,  51,  and  52).  The  global  concept  provides 
a quantitative  spatial  evaluation  rather  than  one  limited  to  the  very 
near  vicinity  of  a point  and  the  subsequent  assumptions  of  homogeneity 
and  isotropy  in  the  case  of  the  classical  theories  (References  11  and 
12). 


3.1  Hypotheses 

The  principal  hypotheses  embodied  in  this  research  are: 

1.  The  response  of  a point  on  the  surface  of  a pavement- 
subgrade  system  when  subjected  to  vehicular  loadings 
can  be  characterized  by  a (Kelvin)  lumped  parameter 
model.  The  three  parameters  of  the  model  are  m,  c, 
and  k,  as  shown  in  Figure  1.  The  equation  describing 
the  hypothesized  lumped  parameter  system  is 

I 

i 

mz(t)  + cz(t)  + kz(t)  ■ F(t)  (1) 
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Inertial  components  of  the  pavement-subgrade  system) , 
c is  the  damping  (the  time  dependent  viscoelastic  or 
energy  dissipating  component),  and  k is  the  stiffness 
(elastic  component);  z(t),  z(t)  and  z(t)  are  the 
equivalent  acceleration,  velocity  and  deflection, 
respectively.  F(t)  is  the  equivalent  input  force  due 
to  the  wheel  load,  (its  maximum  value  for  the  F-4 
aircraft,  for  example,  is  approximately  25  kips). 
Figure  1 shows  the  hypothesized  model. 

The  salient  characteristics  of  the  time  dependent 
response  of  the  pavement -sub grade  system  to  a complete 
vehicular  pass  can  be  obtained  in  a nondestructive 
manner. 

Equivalent  parameters  describing  the  pavement-subgrade 
system  stiffness  (k),  energy  dissipation  (c)  and 
Inertial  (m)  characteristics  can  be  extracted  from  the 
measured  response. 

The  extracted  parameters  accomplish  the  following: 

(a)  they  reflect  the  global  structural 
state  of  the  pavement -subgrade 
system  and;  whereas,  they  are  responsive 
to  seasonal  and  local  variations,  they 
can  be  used  as  the  basis  of  a 
nondestructive  evaluation  scheme. 

13 


•SS, 

m 


(b)  they  can  be  used  to  predict  the 
deflection  response  to  varying 
magnitudes  and  configurations  of 
wheel  loads. 

5.  Many  current  parameters  used  in  conventional  design 
and  analysis  can  be  obtained  from  the  noted  response 
of  a pavement-subgrade  system  to  actual  vehicular 
loading. 

The  developed  theory  and  methodology  provide  the  following:  (1) 
the  theoretical  basis  of  a pavement  evaluation  scheme.  This  includes 
the  determination  of  the  deflection  response  of  a pavement-subgrade 
system  to  prototype  loading,  and  the  methodology  of  extracting 
parameters  describing  the  equivalent  mass,  stiffness  and  energy 
dissipation  characteristics  of  the  pavement-sub grade  system.  (2)  The 
methodology  for  predicting  the  deflection  response  of  flexible  pavements 
to  prototype  loadings  of  varying  magnitudes  and  configurations. 


3.2  Evaluation 

The  pavement  evaluation  scheme  consists  of  determining  the 
parameters  inherent  in  the  transfer  function  of  the  pavement-subgrade 
system.  The  transfer  function  is  a mathematical  expression  which 
represents  the  mechanism  that  converts  the  equivalent  input  to  an 
output.  Figure  2 shows  an  overview  of  the  input-transfer  function- 
output  interrelation  as  it  applies  to  wheel  load-pavement-subgrade 
interaction. 

The  basic  elements  of  the  evaluation  scheme  are: 


(a)  the  signature  determination,  and 
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(b)  the  parameter  determination 


3.2.1  Signature  Determination 

The  signature  is  the  pavement  deflection-time  plot  corresponding 
to  the  pavement  deflection  at  the  outer  edge  of  the  tire  print  as 
shown  in  Figure  3.  Because  the  mechanical  setup  of  the  measuring 
equipment  precludes  the  measurement  of  deflections  closer  than  3 inches 
from  the  outer  edge  of  the  tire  print,  the  signature  is  calculated  from 
deflections  measured  in  the  field  using  the  "stationary  beam-moving 
load"  procedure. 

In  this  procedure  (Figures  3 and  4)  a cantilever  beam  with  its 
pivot  outside  the  deflection  basin  supports  sensing  elements  within  the 
deflection  basin.  The  sensing  elements  measure  relative  motion  between 
the  pavement  surface  and  the  beam  as  the  load  vehicle  travels  in  a 
direction  perpendicular  to  the  beam  position.  The  sensing  elements  are 
positioned  at  different  lateral  distances  from  the  load  wheel  and 
therefore  record  the  deflection  time  history  for  one  complete  pass  of 
the  load  vehicle  at  the  points  of  measurement.  The  location  of  the 
tire  print  is  determined  from  the  Impression  the  tire  makes  on  silver 
adhesive  tape  placed  on  the  pavement  in  its  path  in  front  of  the  beam. 
The  distance  from  the  outer  edge  of  the  tire  print  to  the  first  gage 
is  measured.  From  this  initial  measurement  and  the  fixed  intergage 
distances,  the  distance  of  each  gage  from  the  outer  edge  of  the  tire 
print  is  calculated  [Figure  3(b)].  Figure  4 shows  an  Air  Force  P-2 
fire  truck  about  to  make  a pass  by  two  stationary  measurement  systems. 

The  method  of  calculating  the  signature  is  as  follows: 

15 


Figure  2.  An  Overview  of  the  Transfer  Function  Concept  as 
Applied  to  Wheel-Pavement  Subgrade  Interaction 
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Load  Wheel 


(a)  Profile 


(b)  Plan  View 


Figure  3.  Location  of  Signature  and  Relative  Gage  Location  and  Setup 
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(a)  The  maximum  value  of  the  deflection,  z 


at  each 


! 


particular  gage  (j)  is  obtained  from  each  of  the 


deflection-time  plots  recorded  at  said  gage  (j) 


during  one  pass  of  the  vehicle.  V and  W of  Figure  5 


show  z and  z (z  , z and  z 

raax^  tnaXj  max^  max^  max^  are  not 


shown),  in  perspective. 


(b)  The  distances  x^  [Figure  3(b)]  from  the  outer  edge  of 


the  tire  print  to  each  gage  as  recorded  in  the  field, 


and  the  corresponding  maximum  gage  readings,  z^  , are 


then  fitted  with  a curve  given  by  the  equation 


z ■ A , exp  (B  . x.  peflk) 
maXj  peak  * peak  j 


where  A B^^  a°d  rpeak  are  Parameters  describing 


the  attenuation  of  the  maximum  values  of  the  deflection- 


time plot 8 recorded  on  each  gage.  Equation  (2)  is 


similar  in  form  to  those  used  by  Boyer  (Reference  11) 


and  Hlghter  (Reference  12),  but  has  three  parameters. 


A , , B , and  r , , Instead  of  two  that  they  used, 
peak  peak  peak 


Baladi  (Reference  52)  also  uses  a three-parameter 


equation  to  describe  the  attenuation  of  the  maximum 


lateral  deflection  profile  caused  by  vehicular  loads 


on  highways. 


The  curve  fitting  calculations  are  performed  by  the  computer 


subroutine  EXACFIT  listed  in  Appendix  C.  The  equations  used  are 


derived  as  follows: 
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Given  the  general  equation  of  the  curve  in  the  form 


z m t exp  (0xr) 


(2a) 


This  equation  has  three  unknowns.  A,  6 and  r.  Using  this  basic  equation 

and  values  of  (x  , z ),  (x-,  z ) and  (x-,  z ),  (where  the 

i raax^  z i max^ 

numbers  in  the  subscripts  refer  to  gages  numbered  from  1,  the  closest 
to  the  wheel,  to  5 the  farthermost  from  the  wheel),  the  following  are 
had: 


- AJe.k  ^eak  *1  P“k>  <*■> 


e*P  <6'  *,  Peak)  (3b) 


max-  peak 


peak  2 


‘maxj  " ^peak  exp  <^eak  *3 
Dividing  equation  (3a)  by  equation  (3b)  gives 


‘peak 


) (3c) 


z /z 
max^  max2 


<Vk  ^rp“k  - %..k  »2rpeak) 


Taking  the  natural  logarithm  of  both  sides  produces 


£n  z - in  z 
‘l 


max,  max2  peak 


■ B;.ak  <’l  - *2rpeak>  <«> 


Similarly  from  equation  (3a)  and  (3c) 

Jtn  z - t„  a - 6’  (x.peak  - *,rpeak)  (5) 

max,  max-  peak  1 3 ' ' ' 


Dividing  equation  (4)  by  equation  (5)  and  rearranging  terms  gives 


‘SEEP 


in  z - in  z 


i | Ml  £»  Ml  6 || 

r*  , r*  , f max,  maxo  r*  « r , 

pe*  . x,  »-  - x2  >“k  (6.) 
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Equation  (6a)  is  solved  by  an  Iteration-interpolation  procedure  detailed 


In  Appendix  C,  and  r^^  Is  found. 


^peak  i8  t*len  ca^cu^ate<^  using  equation  (5) , where 


in  z - in  z 

max^  max3 

P6ak  rpeak  _ rpeak 
X1  x3 


Then  A’^^  Is  calculated  using  equation  (3a)  where 


A’  . ■ z /exp  (B1  , x.  Peak) 

peak  raax^  r peak  1 


Usln8  <V  (x3-  ■„>,>  “,1  (V  zmax5>  lnsteod  of  <V  ‘mXj’1 


(x2,  zoax  ) and  (x3,  zmax  ),  equations  (3a),  (3b)  and  (3c)  are  again 


solved  to  find  r"  , , B"  , and  A"  . . The  desired  r . , 3 , and 

peak’  peak  peak  peak’  peak 


^peak  are  ^oun<*  using  the  equations 


rpeak  (rpeak  + rpeak)^ 


6peak  “ (6pcak  + 6peak)/2 


V*  ■ (AP«k + WJ'2 


rpeak’  ^peak  Apeak  are  tlie  Parameters  which  describe  the  maximum 


lateral  deflection  profile,  partly  shown  in  Figure  3 as  the  curve  V'VW. 
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Figure  5.  Overview  of  Method  of  Calculation  of  Signature  Due  to  Prototype  Loading 


(c)  The  signature  is  calculated  using  the  r 


. and  6 . 

peak  peak 

parameters  and  the  basic  equation,  equation  (2a),  In  the  form 


A(t)  - ^(O/exp  (Bpeak  x peak) 


(10) 


where  z^(t)  refers  to  the  deflection  at  gage  1 (closest  to  the  wheel), 
and  t is  the  time  shown  discretized  as  i in  Figure  5.  Subroutine 
SIGNTR  in  Appendix  C performs  these  computations. 

The  signature  thus  calculated  is  the  pavement  deflection  due  to 
prototype  loading  at  point  S in  Figure  3.  Using  this  signature  as 
the  output  (Figure  2),  parameters  k,  c,  and  ro  describing  the  transfer 
mechanism  of  the  three-parameter  model  of  the  pavement-subgrade  system , 
are  calculated. 


3.2.2  Parameter  Determination 

The  objective  of  the  theoretical  procedure  is  to  determine  the 
equivalent  m,  c and  k parameters  of  a pavement-subgrade  from  the 
signature  obtained  from  actual  field  loadings  using  transfer  function 
theory.  These  parameters  define  the  transfer  mechanism  that  converts 
wheel  loadings  to  pavement  deflections. 

The  wheel  load-pavement-subgrade  interaction,  as  shown  in 
Figure  2,  can  be  represented  mathematically  by  the  equation  (1): 

mz(t)  + cz(t)  + kz(t)  ■ F(t)  (1) 

It  can  also  be  represented  by  the  equation 

I(t)  * G(t)  - ^ * G(t)  - z(t)  (11) 


where  F(t)  is  the  equivalent  loading  (input)  function,  m is  the 


equivalent  mass,  G(t)  is  the  reduced  transfer  function,  z(t)  is  the 


response  and  the  symbol  (*)  represents  the  mathematical  operation  of 


convolution  (Reference  53) . In  its  integral  form  equation  (11)  becomes 


IlLl  . 


G(t-T)  dT 


G(t) , the  reduced  transfer  function  in  equations  (11)  and  (12)  is 


obtained  from  equation  (1)  as  follows: 


Consider  the  basic  equation  of  the  model: 


mz(t)  + cz(t)  + kz(t)  - F(t) 


taking  the  Laplace  transform  of  both  sides,  imposing  initial 


conditions  (pavement  is  assumed  to  be  at  rest  prior  to  loading)  of 


z(0)  - z(0)  - 0,  and  transposing  terms,  yields: 


= (a)  * ^ * ( 2 " "k) 

'S  + — 8 + — ' 


where  the  bars  denote  the  transformed  function.  The  reduced  transfer 


function  is  defined  as 


G(s)  - 


2 c k 
s + — s + — 
m m 


Inverting  G(s)  gives. 
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(13a) 


if  k/m  < — =-,  namely,  the  system  is  overdamped , then 
4m 
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Two  ways  to  obtain  the  equivalent  input  (loading  function) 

F(t)  will  be  pursued  in  this  work.  The  first  way  is  by  the  direct 

substitution  in  the  equation 


mz(t)  + cz(t)  + kz(t)  * F(t) 


(1) 


of  the  calculated  values  of  z(t),  z(t)  and  z(t)  obtained  from  the 
measured  deflections,  and  the  parameters  of  the  transfer  mechanism 
m,  c and  k.  The  second  method  is  by  implicit  convolution  using 
equations  (11)  and  (12).  The  implicit  discrete  form  of  equation  (12) 
is  given  (Reference  53)  by  the  equation 


KD  " 


k-i-1 

0(i)/AT  - l I(m) *G(k  + 2 - m) 
m»l 

KD 


(14) 


where  G()  ■ the  reduced  transfer  function 
0(i)  ■ the  output  z(i) 

mx  I(i)  “ F(i)  ■ the  equivalent  loading  function 
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(15) 


AT  - time  increment 


r«.*  -wm  vwv  rm  T«m 


rrw^CTTr^  aJT*n  m jbu"«  im 


and,  the  Initial  value  of  I(i)  ■ 1(1)  is  given  by  the  equation 


1(1) 


. . ,0(1? 

G(l)  «AT 


1 

I 


The  same  two  arguments  above  lead  to  two  ways  of  finding  the  output 
z(t),  given  the  forcing  function  and  the  parameters  m,  c,  k.  The 
first  is  by  the  closed  form  solution  of  equation  (1) . The  second  is 
by  explicit  convolution  using  the  equation 

i 

0(i)  - z (1) m l G(i) * 1 (i  - j + 1)  AT  (16) 

J-l 

where  z(),  G(),  1()  and  AT  are  as  defined  in  equation  (14). 

During  field  testing,  only  the  deflections  with  time  are  measured, 
and  from  these  deflections  the  signature  A(t)  in  equation  (10)  can  be 
obtained  as  discussed  previously.  How  then  are  the  parameters  of 
transfer  frunction  and  the  loading  functions  determined? 

In  concept,  the  methodology  used  is  as  follows:  The  signature 
A(t)  is  taken  to  be  equivalent  to  the  response  z(t)  in  equations  (11) 
to  (16)  as  discussed  in  the  previous  section  entitled  "Signature 
Determination".  The  input  function  F(t)  is  now  known  initially.  It 
is  assumed  however  that  the  maximum  value  of  F(t)  is  the  wheel  load, 
or  approximately  25  kips  (Reference  54)  for  the  F-4  aircraft,  which 
is  the  main  vehicle  of  Interest  in  the  present  study.  Using  the 
signature  z(t),  the  maximum  assumed  value  (25  kips)  of  the  equivalent 
input,  and  the  equations  developed  previously  in  this  section,  the  k, 
c,  and  ra  parameters  of  the  pavement-subgrade  system  are  determined. 
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The  detailed  procedure  for  finding  k,  c,  and  m is  as  follows: 

(1)  The  signature  is  differentiated  numerically  to  yield 
the  time  history  of  the  velocity  z(t)  and  the 
acceleration  z(t).  Typical  plots  of  z(t),  z(t), 
and  z(t)  are  shown  in  Figure  6. 

(2)  The  components  of  the  equivalent  force,  due  to  each 
of  the  three  parameters,  on  the  left  hand  side  of  the 
basic  equation  (1),  are  then  considered  for  three 
special  conditions,  labelled  as  sections  (1),  (2), 
and  (3)  in  Figure  7. 

At  section  (1),  the  point  of  inflection  of  the  signature,  the 
following  equations  hold: 

2(t1)  - 0 
mz(t^)  * 0 

hence, 

cz(t1)  + k zO^)  - F(tx)  (17a) 

At  section  (3),  the  maximum  deflection  response,  the  following 
equations  hold: 

*(t3)  - 0 
cz(t3)  - 0 

hence , 
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Figure  7.  Determination  of  Parameters  of  Transfer  Function  - A Graphical  Overview  of  Methodology 


(17) 


mi’(t3)  + k*(t3)  - P(t3) 


The  third  special  condition  occurs.  Cor  each  deflection  data  set. 


at  a time  between  the  time  corresponding  to  the  point  of  Inflection 


[section  (1)]  and  the  maximum  deflection  [section  (3)].  With  reference 


to  Figure  7,  with  increasing  times  from  section  (1)  to  section  (3), 


the  stiffness  component  [kz(t)l  of  the  equivalent  input  force  Increases 


from  B*  to  a maximum  at  P'.  During  this  passage,  the  damping 


component  [cz(t)]  decreases  from  a maximum  of  A to  zero  (P) . Similarly, 


the  Inertial  component  [mz(t)l  decreases  from  zero  (B)  to  a negative 


value  of  R.  The  time  corresponding  to  the  conditions. 


cz(t2)  - - m*(t2) 


and,  hence,  kz(t2>  " F(t2) 


will  be  called  the  crossing  point.  It  provides  the  final  special 


condition  for  determining  the  measure  of  the  parameters  m,  c and  k. 


This  condition  Is  represented  as  section  (2)  of  Figure  7.  It  should 


be  noted  that  between  section  (1)  and  the  crossing  point. 


kz(t)  < F(t2) 


and  between  the  crossing  point  and  section  (3) 


kz(t)  > F(t2) 


In  concept,  therefore,  the  crossing  point  represents  the  time  at  which 


the  Winkler  hypothesis  (Reference  SS)  Is  satisfied. 


tr 


v 


(3)  Using  the  equation  In  paragraph  (2)  and  the  developed 
method,  an  Iterative  procedure,  given  In  detail  In 
Appendix  C was  used  to  determine  the  k,  c and  m 
parameters  In  this  study.  Briefly  the  procedure  is 
as  follows: 

The  signature  of  the  load  vehicle  is  obtained  as  described 
in  the  section  entitled  "Signature  Determination".  A first 
trial  value  of  k is  taken  as  a fraction  (AA1  < 1.0)^  of  the 
wheel  load  (25  kips)  divided  by  the  maximum  value  of  the 
signature.  Then  it  is  assumed  that  another  fraction  (AA2  < 1) 
of  the  wheel  load  acts  at  the  time  corresponding  to  section  (3) 
Figure  7.  The  flr^t  trial  value  of  m is  calculated  for  the 
conditions  at  section  (3)  using  equation  (17).  The  initial  trial 
value  of  the  c parameters  is  then  computed  using  equation  (18)  at 
the  first  estimated  crossing  point  in  the  data  set,  l.e.,  at 
the  first  discrete  Interval  of  time  after  section  (1).  Using 
these  first  estimates  of  k,  c and  m the  transfer  function  G(t) 
[equation  (13)]  is  had  [EXACT  l]^.  Next  the  loading  function 
F(t)  is  obtained  by  Implicit  convolution  [equations  (14)  and  (15) 
and  NWCONV ] in  the  time  domain  using  the  signature  deflection 
z(t)  is  output  and  G(t).  The  maximum  value  of  F(t)  is  then 
found.  A check  is  made  to  see  if  this  (maximum)  value  of  F(t) 
satisfies  the  necessary  (500  pound)  criterion  that. 


capital  letters  in  parentheses  refer  to  computer  programs  or  symbols 
used  in  computer  programs  in  Appendix  C. 
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24,500  lb  F(t)  maximum  < 25,500  lb 


(20) 


R 


1 

N 

If  this  criterion  is  not  satisfied,  the  next  crossing  point  is  selected, 
c is  recalculated  and  the  same  procedure  is  followed  to  obtain  the 
maximum  value  of  F(t),  which  is  again  examined  with  respect  to  the 
previously  stated  criterion.  This  process  is  repeated  until  the 
possible  crossing  points  are  exhausted  i.e.  section  (3)  Figure  7 is 
reached,  or  the  maximum  value  F(t)  satisfies  the  previously  stated 
criterion,  whichever  occurs  first. 

If  all  possible  crossing  point  times  are  exhausted  a new  value  of 
m is  calculated  by  incrementing  the  value  of  AA2.  Using  the  first  trial 
value  of  k and  the  second  trial  value  of  m a new  c is  calculated  for  the 
first  estimated  crossing  point.  The  procedure  is  repeated,  checking  the 
maximum  value  of  F(t)  against  the  500  pound  criterion.  If  this 
criterion  is  not  satisfied  the  sequence  of  crossing  points  is  again 
stepped  through.  If  these  are  exhausted  a next  trial  value  of  m is 
assumed  and  the  entire  process  is  repeated.  If  the  incremented  value 
of  AA2  reaches  unity,  a second  trial  value  of  k is  calculated  by 
incrementing  Ml,  and  the  entire  procedure,  including  a new  sequence 
of  m and  c values,  is  examined  at  all  possible  crossing  points.  This 
process  is  repeated  until  the  500  pound  criterion  (equation  20)  is 
satisfied. 

If  the  500  pound  criterion  is  satisfied,  a new  c value  is 
calculated  for  conditions  corresponding  to  the  time  at  the  inflection 
point  of  the  signature  (section  (1),  Figure  7,  equation  17a);  using  the 
current  k value  and  the  force  obtained  by  convolution,  corresponding  to 
the  time  of  inflection.  If  this  new  c value  does  not  correspond  to  the 
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c value  which  satisfied  conditions  at  the  crossing  point  and  the  500 
pound  criterion,  the  iteration  process  is  continued  by  selecting 
another  crossing  point  and  repeating  the  entire  procedure.  If  the  c 
values  do  correspond  then  the  corresponding  k,  c,  and  m values  are  the 
desired  parameters. 


3.3  Prediction  of  Pavement  Response  (Signature) 

The  theory  and  methodology  for  signature  and  parameter 
determination  developed  in  the  previous  sections  were  used  to  predict 
the  signature  of  a vehicle  (vehicle  A)  at  a site  prior  to  the  actual 
passage  of  the  vehicle  at  that  location.  This  was  done  by  first 
obtaining  signatures  of  both  a standard  vehicle  and  vehicle  A at 
another  site,  and  the  signature  of  the  standard  vehidle  at  the  new  site. 
This  is  accomplished  by  using  an  equivalency  function. 

The  equivalency  function  is  defined  as  the  ratio  of  the  loading 
(input)  function  of  a vehicle,  F^t),  and  the  loading  function  of  the 
standard  vehicle,  F2(t),  both  obtained  at  the  same  location  (at  a 
standard  site).  This  relationship  is  given  by  the  equation. 


Equivalency  Function  ■ EF(t) 


Fx(t) 

VO 


(21) 


The  procedure  for  determining  the  equivalency  function  is  as 
follows:  with  reference  to  Figure  8,  Block  I,  at  a standard  site  the 
signature  due  to  a standard  vehicle  is  obtained  as  discussed  in  the 
section  entitled  "Signature  Determination".  Using  this  signature  as 
output,  the  transfer  function,  G2(t),  of  the  standard  site  and  the 
loading  F2(t)  of  the  standard  vehicle  at  the  standard  site  are 
determined  as  previously  discussed  in  the  section  entitled  "Parameter 


Determination".  Next  the  signature  of  A (the  vehicle  whose  response  Is 
to  be  predicted)  Is  obtained  at  the  standard  site.  Using  the  transfer 
function  of  the  standard  site  as  determined  by  the  standard  vehicle  and 
the  signature  of  vehicle  A at  the  standard  site,  the  loading  function 
F^(t)  of  vehicle  A is  obtained  by  implicit  convolution  using  equation  14 
and  15,  and  subroutine  NWCONV  (Appendix  C).  Once  the  equivalency 
function  is  had  from  equation  21,  the  predicted  signature  is  obtained  as 
follows : 

First,  the  signature  of  the  standard  vehicle,  z(t),  its  loading 
function  FgA(t)  and  the  transfer  function  GA(t)  are  obtained  at  the  new 
site  as  discussed  in  the  preceding  paragraph.  These  functions  are  shown 
schematically  in  Block  III  Figure  9.  Next  the  loading  function  F^(t) , 
of  vehicle  A at  the  new  site  is  calculated  using  the  equation 

EF(t)  x FsA(t)  - Fu(t)  (22) 

Equation  22  is  represented  in  Figure  9 as  the  Connection  between 
the  loading  function  of  Block  III  and  the  equivalency  function  at  the 
point  labelled  $. 

The  predicted  signature  of  vehicle  A at  the  new  site  is  computed  by 
the  explicit  convolution  of  the  loading  function  of  vehicle  A at  the 
new  site  and  the  transfer  function  of  the  new  site  as  determined  by  the 
standard  vehicle,  using  equation  (16), 

i 

z(i)  - l G(i)  • I (i  - J + 1)  AT 

J-l 

where  G(i)  ■ transfer  function  of  the  new  site 

I(  ) - the  equivalent  loading  function  of  vehicle  A at  the  new 
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site  (see  equations  11  and  15) 


z(i)  m the  predicted  signature 
AT  ■ the  time  increment 
j ■ dummy  variable,  and 

i - dummy  variable  Indicating  the  time  interval  number 
The  explicit  convolution  is  represented  in  Figure  9,  Block  IV  as 
the  connection  between  the  equivalent  loading  and  the  transfer  function 
at  the  point  labelled  y of  Figure  9. 

The  computer  program  PREDCT  that  performs  the  noted  calculations  to 
generate  the  predicted  signature  is  given  in  Appendix  C. 

The  evaluation  and  prediction  methodologies  are  clearly  dependent 
on  the  availability  of  reliable  equipment  for  measuring  pavement 
deflections. 


SECTION  IV 


NEW  INSTRUMENTATION  FOR  DEFLECTION  MEASUREMENT 
Two  mobile  systems  were  developed  during  this  study  for  the 
measurement  of  pavement  deflections:  the  Light  Emitting  Diode  (LED) 

3 

system  and  the  Linear  Variable  Differential  Transformer  (LVDT)  system. 

The  LED  system  is  essentially  an  optical  system  which  does  not  make 
contact  within  the  deflection  basin  of  the  pavement.  On  the  other  hand 
the  sensing  elements  of  the  LVDT  system  make  contact  with  the  pavement 
surface.  A summary  of  the  research  effort  expended  in  the  development 
of  these  systems  and  their  basic  developments  is  presented  below. 


4.1  The  LED  System 

The  LED  measurement  system  consists  of  three  major  operational 
sections:  they  are, 

a.  the  beam 

b.  the  sensor  modules,  and 

c.  the  recording  system 
a.  The  Beam 

The  beam  is  made  of  aluminum.  Its  principal  function  is  to  support 
the  sensors  above  the  pavement  surface.  It  is  designed  so  that  a 
20-mile-per-hour  gust  of  wind  will  deflect  the  end  of  the  beam  less  than 
0.004  in.  Figure  10  shows  schematic  details  of  the  beam.  Figure  11 
shows  the  LED  beam  setup.  Seven  LED  sensor  modules  are  shown  between 
the  channel  of  the  beam. 


3 The  LVDT  system  was  developed  by  Baladi  (Reference  52)  as  part  of  his 
research  effort,  sponsored  by  the  Federal  Highway  Administration. 
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Figure  10.  Light  Emitting  Diode  (LED)  Beam  Schematic  Details 


b.  Sensor  Modules 


Several  ways  of  obtaining  deflections  using  optical,  non-contact 
techniques  were  considered  prior  to  selection  of  the  final  sensor. 

These  methods  can  be  categorized  under  two  general  classes;  namely, 
point  methods  and  global  methods.  The  point  methods  considered  were 
focal  point  displacement,  beam  triangulation,  retro-reflector  beam 
displacement  methods,  diffraction  gages,  interferometers.  Moire  gages 
and  grating  gages.  The  global  techniques  considered  were  holography 
and  Moire  contour  techniques. 

The  technique  selected  for  development  was  the  beam  tri angulation 
method,  using  a light  emitting  diode  as  the  coherent  light  source. 

Figure  12  shows  the  basic  geometrical  optics  of  the  beam  triangulation 
method.  Basically,  pavement  deflection  is  measured  bv  the  amount  of 
displacement  of  the  reflected  light  impinging  on  the  collector  lens. 
Research  has  established  that  a line  of  light  impinging  on  the  surface 
of  small  chips  of  various  materials,  asphalt,  concrete,  and  sand  paper, 
produces  a reflected  signal  that  can  be  normalized  to  produce  a linear 
position  measurement  independent  of  surface  reflectivity.  The  sensors 
have  a range  of  0.2  in  and  an  accuracy  of  0.0005  in.  Figure  13  shows  a 
complete  sensor  module  ready  for  mounting  on  the  beam.  Additional 
details  are  given  in  Appendix  A. 
c.  The  Recording  System 

Signals  from  the  sensors  are  collected  on  a magnetic  tape  unit. 

4 

Figure  14  shows  the  recorder  used  during  this  investigation.  Figure  15 


^A  light  beam  recorder  was  also  used  in  some  special  tests. 
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demonstrates  the  basic  arrangement  of  the  data  on  the  tape.  Additional 
details  on  the  tape  decoding  and  data  reduction  computer  programs  are 
presented  in  Appendix  C. 

It  should  be  noted  that  in  the  LED  system: 

1.  The  output,  after  being  recorded  on  the  magnetic  tape, 
must  be  fed  into  a digital  computer  for  data  reduction 
and  analysis,  and 

2.  No  contact  (excepting  infrared  light)  is  made  between 
the  sensing  elements  and  the  pavement  surface.  The 
LVDT  system  required  some  (minor)  contact  within  the 
deflection  basin  of  the  pavement. 


c.  the  recording  system 

a.  The  Beam 

The  beam  is  a modified  wide  flange  aluminum  I-beam.  Its 
dimensions  and  principal  features  are  shown  schematically  in  Figure  16. 
A picture  of  the  beam,  set  up  in  position  for  measurements,  is  shown  in 
Figure  17. 

b.  The  Sensors 

The  sensors  are  LVDT's.  These  are  fixed  on  the  beam  as  shown  in 
Figure  17,  with  the  sensing  elements  in  contact,  through  adjustable 
jacks,  with  the  pavement  surface.  The  relative  movement  between  the 
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core  of  the  sensing  element  and  the  induction  coil  of  the  LVDT,  fixed 
on  the  beam,  produces  a signal  which  la  a measure  of  the  displacement, 
c.  The  Recording  System 

The  LVDT  recording  system  is  a six  channelled  light  beam  recorder. 
Figure  18  shows  the  recorder  used  in  this  study.  An  analogue  signal  is 
recorded  on  light  sensitive  paper  tape.  The  paper  tape  output  can  be 
converted  to  digital  data  by  using  a digitizing  system  (Reference  57) 
which  automatically  provides  punch  data  cards. 

It  should  be  noted  that,  although  the  LVDT  system  reouires  contact 
between  the  pavement  and  the  sensor,  its  simplicity,  directness  of  data 
collection  and  ease  of  operation  make  it  very  adaptable  for  general  use. 

The  LVDT  and  the  LED  systems  were  both  used  at  various  stages  of 
this  study. 
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SECTION  V 

FT ELD  INVESTIGATIONS 

The  results  of  two  series  of  field  investigations  were  used  in  the 
validation  of  the  evaluation  and  prediction  schemes  discussed 
previously.  The  first  series  of  field  investigations,  Series  A,  was 
conducted  during  the  course  of  this  project.  The  second  series. 

Series  B,  had  been  conducted  in  1972  at  Kirtland  AFB  and  the  basic  data 
were  reported  by  Boyer  (Reference  11). 

5.1  Series  A Field  Investigations 

The  principal  objectives  of  this  series  of  field  investigations 

were: 


1.  To  verify  that  the  new  nondestructive  testing  equipment, 
i.e.,  LED  beam  (and  LVDT  beam),  can  provide  measurements 
of  pavement  deflection  consistent  with  destructive  test 
measurements. 

2.  To  use  the  NDT  (nondestructive  testing)  equipment  to 
obtain  deflection  measurements  for  several  different 
vehicles  at  different  sites,  for  the  purposes  of  providing 
initial  input  into  the  evaluation  scheme  developed  in 
conjunction  with  the  NDT  equipment. 

3.  To  show  that  the  evaluation  and  prediction  schemes  are 
consistent  with  the  factual  information  currently  available. 

The  methodology  evolved  over  the  period  of  this  study  therefore 
reflects  the  need  to  satisfy  the  principal  objectives. 

The  basic  methodology  with  respect  to  equipment  development,  shown 
in  Table  2,  evolved  over  the  period  of  this  study.  As  can  be  seen  from 
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TABLE  2.  FIELD  TESTING  METHODOLOGY  OF  DEVELOPMENT 


Method,  of 

Measuring 

Deflections 

Features 

of 

Method 

Distribution 

of  1 

Measurements 

Gages  Installed 

(a)  Destructive 

(a)  Along  one 

in  Pavement 

(b)  Fixed 

line 

LVDT 

(a)  Nondestructive 

(a)  Along  one 

Beam 

(b)  Contact 

line  and/or 

(c)  Stationary 

Global 

LED 

(a)  Nondestructive 

(a)  Along  one 

Beam 

(b)  Noncontact 

line  and/or 

(c)  Stationary 

Global 

LED 

(a)  Nondestructive 

(a)  Global 

Beam 

(b)  Noncontact 

(c)  Moving 
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Table  2 the  method  of  measurement  of  deflections  progressed  in  four 
steps  from  destructive  fixed  instrumentation,  namely  gages  installed  in 
the  pavement,  to  noncontact,  nondestructive  instrumentation.  Details 
of  the  gage  Installations,  LVDT  beam  and  the  LED  measurement  systems 
are  presented  in  Appendices  A and  B. 

Deflection  measurements  were  obtained  during  each  step  using 
different  vehicles  at  different  sites.  The  developed  instrumentation 
was  used  as  available  for  measurement  purposes  as  shown  in  Figure  19. 
Figure  19  also  shows  a synopsis  of  the  observed  characteristics  of  the 
test  sites,  instrumentation  and  loading  vehicles  used  in  the  Series  A 
field  investigation.  Details  of  the  sites  and  testing  program  conducted 
at  each  site  during  each  phase  are  presented  in  Appendix  B. 

It  should  be  recognized  that  the  sites  selected  reflect  reasonably 
well  the  extremes  of  the  range  of  the  airfield  pavement  conditions. 

Site  1 at  Eglin  AFB  is  a functionally  failed  flexible  pavement  and  Site 
2 is  an  in-service  taxiway.  The  loading  vehicles  used  also  reflect  a 
range  of  tire  pressures  using  airfield  pavements.  The  F-4  aircraft  and 
the  P-2  fire  truck  have  tire  pressures  of  255  psi  and  55  psi, 
respectively.  The  contrast  in  gross  static  wheel  load  also  spans  a 
reasonable  range.  The  F-4  aircraft  and  the  P-2  fire  truck  have  gross 
single  wheel"*  loads  of  approximately  25  kips  and  8 kips,  respectively. 
Figures  20  and  21  show  the  F-4  aircraft  and  the  P-2  fire  truck, 
respectively,  and  their  pertinent'  characteristics.  Figure  22  shows  the 
F-4  load-cart  loaded  and  unloaded.  When  loaded  the  central  load  wheel 
exerts  a load  on  the  pavement  equivalent  to  a single  main  gear  of  the 

^For  the  F-4,  this  refers  to  the  main  gear. 
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EGT.IN  AFB:  SITE  1 


General  Site  Description:  Severely  Cracked: 

used  as  parking  area. 

Date 

Loading 

Instrumentation 

5-75 

(a)  P-2  Fire  truck 

(a)  Gages  installed 

(b)  F-U  Aircraft 

in  pavement 

11-75 

(a)  P-2  Fire  truck 

(a)  Gages  installed 

(b)  F— U Aircraft 

(b)  LVDT  Beam 

(c)  LED  Beama 

3-7  6 

(a)  F-U  Load  Cart 

(a)  LVDT  Beam 

Figure  19.  Synopsis  of  Series  A Field  Investigation 
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(c)  PEASE  AFB 

General  Site  Description:  Uncracked  in-service 

runway:  Rutted  by 
several  passes  of 
KC-135A  Aircraft 
and  P-2  Fire  truck 


Loading 


Instrumentation 


(a)  P-2  Fire  truck 

(b)  KC-135A  Aircraft 


LVDT  Beam 


(d)  EGLIN  AFB:  TAXIWAY 


General  Site  Description:  Cracked  not-in 

service  taxiway 


Loading 


Instrumentation 


F-U  Load  Cart 


(a)  LED  Beam 

(b)  LVDT  Beam 


a Stationary  Beam  - Moving  Load  operation 
b 6000  feet  of  runway  were  also  tested  on  this  date, 
c Moving  Beam  - Moving  Load  operation 

Figure  19.  (Concluded) 
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Gear  Configuration  Single  Wheel 

Test  Weight  25  Kips 


Figure  20.  F-A  Aircraft  and  Pertinent  Characteristics 
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Twin  Tandem 


Twin 


Twin 


STD.  TEST  WEIGHT  66,000  fc 


Figure  21.  P-2  Fire  Truck  and  Pertinent  Characteristics 
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F-4  aircraft.  The  majority  of  testing  was  done  using  this  vehicle. 
Because  of  the  limited  number  of  prime  movers  that  were  made 


available  during  the  periods  of  testing,  it  became  necessary  to 
utilize  some  test  results  previously  reported  by  Boyer  (Reference  11). 

5.2  Series  B Field  Investigations 

The  field  investigations  in  Series  B were  conducted  by  Boyer  at 
Kirtland  AFB  in  New  Mexico  in  1972.  In  his  investigation  deflection 
measurements  were  made  using  gages  installed  in  the  pavement  in  a 
manner  similar  to  the  procedure  used  in  the  Series  A field  investiga- 
tions. For  details  on  his  instrumentation  and  methodology  the  reader 
should  refer  to  the  original  reference  (Reference  11).  The  pavement 
characteristics  at  the  three  sites  reported  are  shown  in  Figure  23. 

The  prime  movers  and  their  characteristics  are  shown  in  Figures  24 
through  26. 

As  previously  stated  the  series  B data  provide  additional 
information  required  to  test  the  evaluation  and  prediction  schemes. 

The  gear  configuration  of  the  C-130,  C-135  and  C-131  provide,  in 
addition  to  a variation  in  gross  load  and  tire  pressures,  a loading 
function  due  to  multiple  and  sequential  wheel  loads,  and  combinations 
thereof.  The  C-131  (single-twin)  is  an  example  of  a multiple  but  not 
sequential  loading  function.  The  C-130  (single-tandem)  is  an  example 
of  a sequential^  but  partially  multiple,  loading  function.  The  C-135 
(twin-tandem)  is  an  example  of  sequential  and  a multiple  loading 

^The  P-2  fire  truck  is  an  example  of  a purely  sequential  loading 
function. 
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Gear  Configuration 
STD.  Test  Weight 


Single  Tandem 
90  Kips  (Gear  load ) 


Figure  24.  C-130  Aircraft  and  Pertinent  Characteristics 
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Gear  Configuration  Twin- Tandem 

STD.  Test  Weight  210  Kips  (Gear  load) 


Figure  26.  0135  Aircraft  and  Pertinent  Characteristics 
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The  data  from  the  Series  A and  the  Series  B Field  investigations 
were  analyzed  using  the  evaluation  and  prediction  schemes  previously 
presented. 
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SECTION  VI 
RESULTS 


This  research  effort  encompasses  a broad  spectrum  of  engineering 
activity  ranging  from  the  development  of  NDT  equipment  and  a procedure 
for  nondestructively  evaluating  flexible  pavements,  to  the  development 
of  a methodology  for  predicting  the  deflection  response  of  various 
pavements  to  different  aircraft.  This  section  presents  the  following: 

1.  typical  results  of  equipment  field  verification  tests, 

2.  results  of  the  evaluation  of  the  two  sites  (Series  A: 

Field  Investigations)  using  the  developed  evaluation 
scheme , and , 

3.  results  of  the  application  of  the  prediction 
methodology  to  three  different  aircrafts  at 

three  different  sites  (Series  B:  Field  Investigations). 

6.1  Equipment 

It  is  necessary  during  the  development  of  new  eouipment  for 
measurement  of  pavement  deflections  to  show  that  the  new  equipment 
provides  measurements  which  are  consistent  with  those  obtained  using 
established  measurement  systems.  It  is  assumed  that  the  I.VDTs 
installed  in  the  pavement  comprise  the  absolute  measurement  system. 

6.1.1  The  LVDT  System 

The  verification  check  of  the  LVDT  system  consisted  of  a 
comparison  between  pavement  deflections  measured  by  an  LVDT  mounted 
on  the  beam  and  an  LVDT  installed  in  the  pavement  as  shown  in 
Appendix  B of  this  report. 

Figure  27  shows  a comparison  of  the  pavement  deflection  measured 
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by  an  LVDT  installed  in  the  pavement  and  an  LVDT  mounted  on  the  beam, 
and  produced  by  the  same  pass  of  a P-2  fire  truck  during  the  testing 
sequence  conducted  at  Eglin  AFB,  at  Site  1 on  November  17,  1975. 

The  results  show  that  the  installed  LVDT  and  the  LVDT  mounted  on 
the  beam  give  basically  the  same  measurements. 

6.1.2  The  LED  System 

The  verification  checking  of  the  LED  system  consisted  of  two 
phases:  a laboratory  and  a field  testing  phase. 

The  laboratory  phase  was  conducted  by  Science  Applications 
Incorporated,  and  is  summarized  in  Appendix  A of  this  report.  The 
principal  field  checkout  consisted  of  a comparison  between  the  output^ 
from  an  LED  on  the  LED  beam  and  the  output  from  an  LVDT  mounted  on  the 
LVDT  beam. 

Figures  28  and  29  show  the  pavement  deflections  measured  by  the 
LED  system  and  the  LVDT  system,  respectively,  and  produced  by  one  pass 
of  a P-2  fire  truck  (Figure  21),  during  the  testing  sequence  at  Eglin 
AFB,  Site  1,  on  November  19,  1975.  Figures  30  and  31  show  the  actual 
digital  output  from  the  LED  system  and  the  analog  output  from  the  LVDT 
system  as  recorded  by  the  light  beam  recorder.  Since  the  deflections 
shown  were  measured  for  two  separate  passes  of  the  P-2  fire  truck,  a 
direct  comparison  of  measured  deflections  can  not  be  made. 

As  stated  previously  a direct  comparison  of  the  measurements  made 
by  the  two  beams  (LED  and  LVDT)  could  not  be  made  since  it  is 

^ Due  to  persistent  problems  of  compatibility  between  the  LED  modules 
and  the  magnetic  tape  recording  unit,  the  LED  output  was  recorded 
on  a light  beam  recorder. 
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Deflection  LVDT  on  Beam  (ft) 


Figure  27.  Comparison  of  LVDT  on  Beam  and 
LVDT  Installed  in  Pavement 
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Figure  30.  Actual  Digital  Output  From  LED  System  On  Light  Beam 
Recorder:  P-2  Fire  Truck  Is  Loading  Vehicle 
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Dots  indicate  points 
used  in  digitizing  light 
beam  output 


Figure  31.  Analog  Output  From  LVDT  System:  P-2  Fire  Truck 


impossible  to  get  the  two  beams  over  the  same  point  at  the  same  time. 
Thus,  a comparison  was  made  of  significant  points  on  the  measured 
deflection  plots;  e.g.,  peaks  and  troughs.  The  points  P,  0,  X and  Y 
were  selected,  where  P is  the  minimum  deflection  between  the  second 
two  peaks  of  the  deflection  time  plot;  0 is  the  minimum  deflection 
between  the  first  two  peaks;  X is  the  first  peak  and  Y is  the  minimum 
between  the  first  two  peaks  and  the  last  two  peaks.  The  ratios  P/Q 
and  X/Y  are  thus  dimensionless,  and  are  Independent  of  the  system 
(dependent)  calibration  constants.  That  is, 


X Actual  Light  Beam  Recorder  Output  at  point  X times  Constant 
Y Actual  Light  Beam  Recorder  Output  at  point  Y times  Constant 


Pavement  Deflection  corresponding  to  X 
Pavement  Deflection  corresponding  to  Y 


Figures  28  and  29  show  that  the  corresponding  ratios  have 
essentially  the  same  values.  Both  the  P/0  and  the  X/Y  ratios  differ  by 
only  0.1. 

From  the  discernable  deflection  observable  at  point  Y in  Figure  28 
and  Figure  29,  it  can  be  seen  that  measurements  of  less  than  0.0001  feet 
are  obtainable  by  both  the  LED  and  LVDT  systems. 

In  June  1976,  both  the  LVDT  and  the  LED  systems  were  used  in 
connection  with  a study  to  investigate  the  difference  in  the  dynamic 
response  of  aircraft. 

A mechanical  system  was  designed  to  allow  the  LED  beam  to  be 
rigidly  mounted  on  the  F-4  load  cart.  The  system  was  fabricated  at 
Tyndall  Air  Force  Base  by  Air  Force  personnel  and  was  then  transported 
to  Eglin  AFB.  At  Eglin,  2000  feet  of  taxiway  and  1 1/4  miles  of  runway 
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were  tested  using  the  LED  system  attached  to  the  F-4  load  cart  as  shown 
in  Figure  32.  The  LED  measured  deflections  due  to  the  load  wheel  of 
the  F-4  load  cart  as  the  load  cart  moved  at  creep  speeds  over  the 
pavement.  LED  measurements  were  made  continuously.  The  LVDT  system 
was  used  to  make  deflection  measurements  at  stations  100  feet  apart  as 
the  load  cart  passed  In  front  of  the  stationary  beam.  The  continuous 
measurements  made  by  the  LED  were  synchronized  with  the  LVDT  stations 
by  the  use  of  an  electronic  event  marker  of  the  LED  system.  The 
synchronization  allowed  the  time  at  which  the  load  wheel  passed  the 
LVDT  stations  to  be  recorded  on  the  magnetic  tape  to  within  0.25  sec 
as  follows.  At  the  instant  the  load  wheel  passed  the  LVDT  station 
marked  on  the  taxlwav,  a button  was  depressed  on  the  control  panel  of 
the  LED  system.  This  caused  a digital  record  to  be  made  on  magnetic 
tape.  Since  the  scan  rate  was  four  scans  per  second,  an  event  marker 
could  be  recorded  to  within  0.25  seconds,  which  constituted  one  scan. 
Figure  33  shows  a summary  of  the  taxiway  deflection  data.  The 
deflections  on  the  dashed  line  in  Figure  33  were  measured  at  locations 
offset  by  about  20  feet  from  the  previous  measurement  locations. 

g 

Listing  of  typical  measurements  made  by  the  LED  system  is  given  in 
Appendix  D. 

The  results  shown  on  Figures  28  and  29  indicate  that  the  pavement 
deflection  measurements  with  the  LED  system  are  consistent  with  those 
of  the  LVDT  system,  in  form;  and  the  LED  can  measure  deflections  with 
the  desired  accuracy  in  a field  environment. 

8 

The  LF.D  system  data  reduction  software,  a complete  listing  of  all 
the  taxiway  data,  and  a copy  of  the  data  (on  magnetic  tape)  has 
previously  been  supplied  to  the  Civil  Engineering  Center,  Tyndall 
AFB,  Florida. 
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6.2  Evaluation 

Evaluation,  as  presented  here,  consists  of  the  following: 


1.  determining  the  parameters  Apeak»  ^peak*  anc*  rpeak 


of  the  maximum  lateral  deflection  profile  due  to 
F-4  loadings, 

2.  obtaining  the  k,  c,  and  m parameters  of  the  transfer 
function,  and 

3.  relating  the  determined  parameters  to  performance. 

The  two  sites  in  Series  A field  investigation  are 
compared  using  the  data  obtained  during  the  field 
testing  program  conducted  at  Eglin  AFB  In  March  1976. 

As  the  F-4  aircraft  is  the  vehicle  of  primary  interest 
in  this  study  only  the  results  obtained  using  the  F-4 
load  cart  are  presented  in  this  section.  A more 
complete  listing  of  the  data  base  and  results  is 
given  In  Appendix  0. 

Briefly,  the  evaluation  procedure  consisted  of  measuring  pavement 
deflections  using  the  LVDT  beam  as  the  F-4  load  cart  passed  in  a 
direction  perpendicular  to  the  stationary  LVDT  beam.  From  these 


deflections,  the  parameters  Bpgg^*  and  r ^ of  the  maximum 


lateral  deflection  profile,  and  the  signature  were  determined  as 
discussed  In  detail  In  paragraph  3.2.1.  The  k,  c,  and  m parameters  at 
each  site  were  then  obtained  using  the  methodology  of  paragraph  3.2.2. 


6.2.1  Pavement  Deflections 


Figure  34  shows  a typical  set  of  deflection  functions  due  to  the 
F-4  load  cart  as  measured  by  the  LVDT  beam.  Gage  1,  which  recorded 
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Figure  34.  System  Response  Functions  F-4  Loading:  Eglin  Site  1 (Typical) 


the  largest  deflection  is  located  closest  to  the  wheel.  Figure  35 

shows  the  maximum  deflection  of  each  gage  plotted  against  the  distance 

of  the  gage  from  the  edge  of  the  tire  print.  The  dashed  line  through 

the  points  G-l  throtigh  G-5  is  the  lateral  profile  of  maximum 

9 

deflections  at  Site  1 for  one  pass  of  the  F-4  load  cart.  Figure  35 

also  shows  the  A , , 3 , . and  r , parameters  determined  as 

peak’  peak’  peak 

discussed  in  paragraph  3.2.1  and  the  maximum  deflections  at  each  gage 

calculated  by  equation  2.  Figure  36  shows  a similar  plot  for  one  pass 

of  the  F-4  load  cart  at  Site  2.  Tables  3 and  4 give  the  stations^ 

and  values  of  the  A , , 3 , . and  r . value  for  Sites  1 and  2, 

peak’  peak  peak  * 

respectively.  Figures  37  and  38  show  the  variability  in  the  measured 
maximum  lateral  deflection  values  for  Sites  1 and  2,  respectively. 
Figures  39  and  40  show  a comparison  of  the  local  and  site  characteris- 


tics at  the  two  sites  as  given  by  the  3peak»  and  r ^ values. 


"Local"  refers  to  a station  or  location  at  a "site".  The  lines  in 


these  figures  are  drawn  for  convenience  and  do  not  represent  continuity 
of  parameter  values.  Figure  41  shows  the  typical  effects  of  two 
consecutive  passes  on  the  maximum  lateral  pavement  deflection  profiles. 
These  results  show  the  following: 

1.  The  gage  (Gage  1,  Figure  34)  closest  to  the  wheel  records 
the  maximum  deflection  and  retains  the  most  permanent  set 
after  the  wheel  load  has  passed.  The  gage  (Gage  4,  Figure 
34)  farthest  from  the  wheel  does  the  opposite;  namely. 


Sites  are  indicated  in  Figures  B-l(a)  and  B-3(a). 

0 

Stations  are  given  in  Appendix  B,  Figures  B-l(b)  and  B-3(b). 
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Figure  35. 


Eglin  AFB:  Site  1 Lateral  Profile  of 

(T^rirflection8  for  f~4  L°ading 
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Figure  36.  Eglin  AFB:  Site  2 Lateral  Profile  of 
Maximum  Deflection  for  F-A  Loading 
(Typical) 


1 


fSIStl&ll 


illllim 


mm  i 


wimimi 


imm  «f 


I IM'MMM 

« * * • « » - a a - 


M 5 MMM 

• * I ****** 


!!!!!!!!! 


I (1MMMM 
1 


!!!!!!!!!! 


!!!!!!!!!! 


1 M3!!?!{  !! 


MMMI  • M 

4*  O M * * * * *• 


s ! n m ! n 


tr  m iiwwri. 
npn  M(k) 


mt 


1 


IU 


to. 


I 1 ! 


Uf 


c • (Ktj)  - 


Civil 


•it*.  la  reality 


(c#, 


• • • 


Y__).  gkmm  bf  tbm 


141 


Nnmm  «•  iMlpai,  kaatecU?,  to  IMt  Mbpraia  MUtttoM. 

«•  WMWtotoi  tm  «Mfe  MbpraiM.  <MU  m tto  Ntov  toto,  tUaait 
Mttltoi  an  Wilt  war  tti«t  atobpraiaa.  Tto  aac  affast  la  dm 


tofUttlMM  W tW  NltlM  «•  ftttoi  U *up  vtUto  tto  MM  Ntot 
•f  aipttato  af  Mil  itvalM  by  tW  iaalpi  pwaaP—a.  tlaaa 


ItftM 


c • MmafM 

> • af  akaalf  af  imc  •(  atmaft 

T • Traffic  «Udi  la  fact  (iwwally,  T.$) 
w • MUtfe  af  cftva  t— tart  araa  af  aaa  tiia  la  tadtaa.  aatf 
I • Cyclaa  af  apavaciaa  (aaa  crda  caaatata  af  aaa  laaiiag  aatf 
aaa  tafcaaff) 


I 

I 


MS 


4AM 

p/< 


Uural  4latrik*c  1m  •(  aircraft  am 

at  U>  Milan  4l«<  rim  in.  t» 

•f  aircraft,  fta  |an  tkt  fallaviaa 


pit  • jr 


U'  la  tto  wi4tk  af  tka  tlra 


par  lack  af  wl4tk,  per  aircraft 


Mick 


>1  mt  aril 

ta  at  tka  palrt 

af  aaataaa  accMalMiaa.  flpr*  Wft)  illaatrMaa  tka  pic  rati* 
aa  apallaf  u a alasla  Maal  ta  aaacapt . witk  raf*rawca  ta 
70(k),  aalp  tka  af  facta  af  tlra  prlati  ar  partial  tin  prlati 
«pplia4  ta  tka  aarfaca  witkia  tka  «ri4tk  U(  ara  caaat4at«4  la  tkia 


Ml  prac  a4ar»  m4  kalaa«  a prepared  aatka4  katk 
tka  tin  prlat  (caw  tact  ara  a vi4tk)  aa  a turp  aaraweter  ia  tka 
4af iaitiaa  af  tka  pic  ratla  Tka  aatka4  4aaaiapai  kara  aaaa  tka 
■Mfa<  fyaawlc  lataral  4aflactiaa  preflle.  4aa  ta  tka  pratatppa 
leaf,  ta  4afiaa  paaa  ta  cmripa . ami  ralMaa  tkia  rat  la  ta  ptrforwrt 
aitaailai  tka  park  eaawlatin  4af lactlaa  raacept  iatraiweaf  kt  Rlphtar 
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V V ...  AV  . - » -V  V - A'J 


(heferewre  12) 


t.2 


IL1- 


Is  f 
1 

CIm  a* 


(aiMsui  cMc#t  «4*rtylac  the  4>flscilwlMii  p/c  raila 
m cIm  follewiag  rtsasslai  The  1ml  sf  ilplflcm*  ef  a 
il  dvr lag  »m  fSM  mr  « pwms t mctlm  is  MlcsM  by 
lateral  tmtsli  4«flact  laa  ksals  Am  the  lateral  dletmes 
free  tike  leas  when.  increases,  the  wart  4am  by  tie  lea4  4ecreaaee ■ 
Heme,  lanadlatel*  jrWe r che  wheel  the  la f lessee  la  greatest ; whereas, 
at  a pel**  where  the  aeaeured  4ef 1 act lea  la  sers , the  laf Isaacs  la 
aetllglble 

As  eas  peas,  la  effect,  produces  4efl set leas  free  a aaalaae 
wa4er  the  wheel  ts  aero  at  sew  41 st sees  free  the  wheal;  far  a p/c 
ratio  to  he  assn lawful , it  shoul4  sc  covet  aot  only  far  the  Ia4aca4 
eastiT  level  corresponding  to  the  as  alma  4sf  lectlea  hot  far  all  the 
levels  la  the  range  The  prlaclpal  raaasalag  hare  is  that  the 
pevamac-  eobgrade  erst  so  rawsahars  1 ast  only  the  aarl—  4ef  lectlea 
hot  all  lea4  la4oce4  deflections . Haaca , la  general.  large  lse4 
Ia4ace4  4efl set  loos  si  s potat  result  la  greeter  oarh 

reeaahraaces  of  It.  thar  4e  son ller  loe4  Ia4aca4  4aflectl 

With  refereace  ts  Figure  71,  each  of  the  corves  label ls4  P,  Q, 

R.  S represents  e series  of  peak  lateral  4ef lect loos  If  see  peas 
at  s laa4  vhee.  Is  applied  at  s pelat  corraapoa41ar  to  % • 0,  the 
aaalwsa  deflection  (at  » • 0)  would  he  A • i#.  If  che  sans  1« 
hewn  applied  at  i • 1,  « asslwaa  deflection  at  1,  would  have 

as  sow red  at  * • 0,  If  the  lead  had  been  applied  at  * • 2.  then 

* 
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Figure  71.  OvtrrlM  ®f  eh*  C**c*pt  U**4  In  Ch*  D®fl®rtlon 


would  um  kM i aiMNi  M i • 0;  similarly  for  a Uai  at  i • ),  see. 

Jm  soeerage  th«,  la  camcagt , la  eta  r— letloe  pMk  daf lact lor 
at  tha  polat  la  quaatloa  chat  earreagmada  to  aulaa  defloct lea  A 
talar  tha  wheel.  far  example,  If  oaa  pees  of  a load  wheel  la  aada  at 
i*0,  thaa  oaa  coverage  reeui’.ts  at  that  palat.  3a  tna  othar  head, 
oaa  pass  aada  at  a • 3 corra agenda  to  tha  fractiaaal  rover age  A 

at  * - 0. 


Suppose  that  f#  through  f q are  tha  dlaerata  freqwamclea  at  i • 0, 
1.  2.  3.  a,  of  tha  dlatrlhatloa  of  passes  of  aircraft  laterally 

shout  * - 0.  Thaa  tha  cumulative  peak  deflect lea  at  i • 0 due  to  thla 
f requeue v distribution  becomes 


at  lee  feak  Reflect loe 


9 lie 


A,  f. 


where  *•  Is  the  nuaher  of  passes  a ad  A,  era  tha  peak  deflect  leas 

p 1 

correspond lot  to  tha  paints  In  quantise.  Tha  coverages  H daa  to 
tho  t passes  la  than, 

-amulet lea  »sak  Oaf lact laa 

c * iaxlua  Oaf  lection  umdar  tha  dhoel 

• % j *i  V4. 

1-u 

The refers , 


q 

JL 

q 

c 


£ 

e 


A 

o 


A1  f 


1 
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Tabl*  7 


tear  41 aerate  41etrihotloee  *r« 
elagle  «kMl  1 aa4.  Figaro  73 
Lateral  41atrlb«tloee  (keya4  by 


la  Figaro  72  far  aa 
the  affect  of  the  fear 
r to  Figaro  72)  aa  the 


ratio 


4af lactlaa.  Figaro  74  above  the  affect  af  the  scatter  of  the 
fregaoeev  41atrlhatlaa,  ae  glvea  by  Ita  ataa4ar4  4aalatlao 
rhaae  raaalta  ahav  that. 

1 the  aeoaapt loo  of  a oalfora  41etrlhatlaa  af 
(41atrlhatlaa(^  ea  Figaro  72)  rlol4a  the  1 
caaalatlve  pask  4af loct loa  gar  paM 

2 aa  weal4  bo  o^acto4 , the  41etrlbetlea  with  the 
aval loat  acattor  elol4a  the  greateet  caaalatlvo 
peak  4eflectl 


taa 


Pf  * 


of  cavalat  lvo  p««t  4of lactlaa 
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m the  Uuttl  iUlrlhtUi  ef  MNM  Mi  tto  pMfc  MUttlM  pnflU. 
the  p/c  rati*  that  pwtiiii  a •(  <anlf?tn  tit  war  Utility  ef 

tit  Ufat  (aircraft  Laaia)  cheractertetlc*  aa*  pevaaaat-eafrgrai* 

ttitli  ra  fa  react  t#  Ptpn  75.  lapet  veriafrlee  caaalat  af  aircraft 
laaia  aa*  lateral  *letrlfretl*a  cherarteriet ice  ef  leaie.  Ttaaa  cea 
iacleie  via*  af facta,  aarface  oaafittea,  type  af  aircraft,  aaf  peaatlty 
ef  faal  afreet*.  Tfra  eat  a act  pet  variable  la  tfra  paafr  pavaaaat 

ef facta  aaf  etfrar  aafrleat  oaeiltiaaa.  Ft  par*  75  efreee  free  tfra  p/c 
rat  la  cm  fre  axillae*  te  eetiaate  tfra  aaeeac  ef  aaargy  Uparte*  te  a 
p*v— at -ivfr  grate  eyetaa  tea  te  variafrle  lapeta  ami  alae  free  te 
eatiaata  the  aefrer  ef  paaaaa  rapilial  te  laieai  a carta  la  aaargy 
lap at  aa  glvea  fry  tfra  tetal  ceaalattve  peafr  *ef lact&eea. 

Far  taaaple,  aaaeaiag  that  the  lateral  flat rlfret lee  af  paaaaa  free 
a etaafari  4edatlea  Of  2.0  (iiatrlfretlea (J) ia  Figara  71),  tfra  aaafrer 
ef  paaaaa  rep  el  re*  te  reach  the  threahel***  tetal  ceaalattve  paafr 
*ef lectlaa  valee  ef  2200  feat  la  1.140,000  paaaaa  a*  500,000  cyclee  ef 
eperatlee. 

Qa  tfra  ether  free*,  if  the  aaeeaa*  lateral  iiatrl fret  lea  af  paaaaa 
free  a ataatfer*  ievletiaa  ef  2.04.  tfraa  500,000  cyclee  af  eperatlee 
aval*  rial*  a ceBaletive  paafr  Oeflectlea  valee  ef  afreet  1500  feat, 
ehlch  la  eall  fre lee  the  tfrreehel***  vale*  far  that  *letrifretlea.  tar 
tfrla  caee  050,000  cycle*  ef  eperatlee  ceel*  fre  telereta*  frefere  aavera 

tsrvehel*  valve  *f  2200  feet  va*  ttv*e  frv  Hight*r  (Oef*reec*  u ) 
vithevt  r**ar*  te  tSe  lateral  pealtiealec  of  traffic,  ^acklag 
*4*ttienal  tnfemattee  at  tai*  tie*  the  tee*  vale*  va*  uee*  t*  the 
pr**eet  ete*v 


154 


I 


4Uim»  atcht  k«  oapoctto 

to  awtktr  •uvU,  Hpf*M  that  4m  to  «*trltr  or  rtMrftdai 
of  MMMtt  , HTMtCtt  a*4  M»  4*(Uttl«M  clua|*  to  tlMM  |1«W  i» 
TakU  I t(Ur  290,000  crcUt  of  •hmUw.  tkw  T 00, 000  cvcloo  a f 
•f«ratt«u  (raclMr  than  290,000)  «m)4  mv  to  roqnlrto  to  *tUU  t to 
thnutll^  ««1«  of  2200  foot.  Tkli  ummi  t tot  t to  Md  1m 
"rwMkan  t to  lotolng  rrtar  to  rtnrfMlii  If  rmrfatlai  tod  mt 
trtw  290,000  tfclti  of  tfaratlo  «mU  tovo  >»w  r*t«l r*4 

to  stow  offsets  of  dttrtM 

Tto  otooo  tutfltf  mw  sttoaptto  to  ihw  ttot  KtMl  — ooorto 
distribution  iota  such  oo  ttot  provldto  by  Motong  'tofarooca  57),  <M 
net no 1 nsoowra*  inflections  which  account  for  eon loot  aoooooal 
offsets)  con  to  root  Into  to  oattnota  tho  conlatlw  toast  a affacta 
4m  to  prototrpo  lotos  to#  a atotoori  alta  ato  atotooi <t  vohtclo 
tovo  toon  aolactto.  tto  iasaaga  off  acta  of  onflows  nognlttoaa  of  lotoa 


StCTlfW  IX 

WMMr 

TW  ftilUji  *1  tlM  yniit  r»M«rth  tfftrt  c—  ho  mmshm 
— 4ar  Ur»«  k«atc  ctc«(«rlM,  mmIt 

I Maitit met l **  M««Ht  uttut  o*«ip— at  4ml*f«nt 
2.  — intractlw  p«w»t  katevtar 

pra4tctlaa,  «i 

ralata4  e*o4*pta. 

T—  afilpHit  afttan  for  aaa  U tka  aaaiaatrattlva  avalaatlaa  af 
flailkla  HMM*ta  haaa  kaaa  4atalaH  > Mr*  tka  LtZ)  arataa  aa4 

tha  LW  ayataa  Th*  (aaaUility  of  ra^Uly  — aaurlag  pa— t 4*f loc- 
tlo—  coo and  kf  aarlac  laaia  oalai  th*  faaalapaf  a^alfoaat  kaa  Woo 
i— ottfto4  ftaU  toot  loo  of  tha  imlapaf  ay  to—  aaui  protKfpo 

Ionia  — o ea«4actof. 

A ooo iaotrocti—  pooo— nt  avalaattoa  aitho4olor  —inf  tra—  for 
faacttaa  thaorr  Wo  kooo  4— lepaf  fara— tor*  of  tho  tra— for 

foaaeti—  ahlck  rofloet  ait*  ckoroctoriat tea  ha—  k—  ioo iato4  mmi 
ralatoi  to  parlor— c*  of  t—  altoa  Para— tor*  of  th*  iatoral  pooh 

ioflocti—  prof  Ho*  ha—  boon  laraotltatai  Slapiifloi  procWor** 
for  *at lot l«t  th*  para— t*r*  of  th*  tra— f*r  faatt too  — 4 para— tor* 
(CM,  ate.)  «aai  to  A** l an  am  kaaa  aauoatoi 

A — thoiolot?  haa  ho—  pro*— 1*4  for  prWictio*  th*  iofloctioo  if 
th*  roap—  at  a ataadari  ait*  t*  a ataaiard  lead  — hlei*  aad  anotnor 
— hicl*.  a—  th*  r**pen—  of  th*  othar  att*a  to  th*  ataaiari  —Kiel* 
ar*  k—  Th*  aa—  tra— f*r  f— eti—  para— t*ra  that  ar*  aaad  ia  th* 


1*0 


sccnr*  x 


cmcutsiom 

This  nmreh  effort  Iim  tttMftM  to  cower  e epee  tree  of 
activities  teat  carries  nondestructive  test  let  aed  mlsttloc  of 
flexible  el r fie  14  psvnntt.  The  following  conclusions  bees  bees 
rosebud  es  e result  of  this  Italy  . 

1.  lap  14  nondestructive  neeeurenaets  of  pavement  4e  fleet  lees 
4ee  to  srrlai  prototype  loe4s  cee  be  obtained  is  tbs  field 
salat  tbs  LO  sts  tee  "brvodboerdsd  is  this  research  effort 
ee4  the  csepealee  LYOT  e vet an 

2.  The  parameters  of  both  the  trees for  fwset  lea  serf  tbs  peak 
lateral  deflection  profile  reflect  elmrec ter lattes  of  tbs 
pass  east -sob grade  spot  one  ae4  can  be  seed  as  tbs  basis  of 
s aoe4ee tractive  paveeeet  evaleutlee  abse.  Sees  csoan 
parameters  such  as  CM  and  nodules  of  sobprads  reaction  can 
be  obtained  free  theee  defined  is  this  study . 

1,  The  developed  eethodolof? . us lag  tbs  sane  paraesters  of 
the  transfer  f see ties  that  arc  seed  to  evaluate  s s evens at 
svotse,  ess  be  seed  to  predict  the  deflection  response 
fextlec  of  s site  to  different  uepnltudos  and  csef  lgwret lei 
of  loadlnas 

4 The  concept  of  cenulst lev  peak  deflection  (Reference  12)  as 
It  relates  to  energy  Induced  Is  « pevunent  subprade  srstsn 
due  to  load  repetition,  was  set ended  to  account  for  tbs 
lateral  distribution  of  aircraft  The  frequency  character- 
istics and  tbs  mader  of  passes  of  aircraft  ere  shewn  to 


SECT!**  X! 


trco'tfni’iATiows 

Thla  wort  w affo rded  the  ^fortualtT  to  'ontlikr  la  mm  depth 
a ti(Rkfi  m portion  of  t!u  car rant  state  of  the  art  of  mstde  a tractive 
tasting  taon  said  net nods  t«M  on  airfield  Mvsarnta  haaed  on  tha 
adtat*  of  thla  sort  the  following  rtcoandit  ions  ara  forward  ad 

1 Tha  (iM  of  prototype  loadings  la  evaluation  should  b« 
encouraged 

2 Tha  lesalijpamt  of  a * 'ast  ( lased  os  t)«  jse  of  nondest  rue 
tivelv  aaaaured  osar  amt  deflections  dws>  to  prototype  loads, 
similar  to  that  brief  Is  disc  us  sac  In  pararraph  4.4  of  this 
wort  should  'ha  pursued 

1 ,v*a  prototype  of  noncontact  measuring  equlpsaent  shown  to  ha 
faaaihla  In  this  Invest i gat  Ion  should  ha  developed  This 
level i twar.t  should  taka  adv an* age  of  the  carrant  atata  of  tha 
art  la  optical  design. 

• A tSortHjg  st odv  of  tha  lateral  distribution  of  aircraft  has 
ret-antly  S**ti  made  i Reference  A similar  lnfomatloa 

rat  taring  st udv  should  be  Initiated  to  eollact  data  an 
•umbers  and  tepee  *>f  aircraft  using  specified  ptvtssstt  at 
bat  i nmer'lai  and  all  tart  .natallat Iona  favenent 
dcflactlona  should  ha  obtained  rarularlv  at  the  sect  *ona 
e ng  monitored  TMs  data  aae  la  nacesaarr  if  tie  concepts 
sid  the  pr*  to  Ine  rt  one 1 us Ions  reached  at  thla  tin*  are  to 
* validated  and  general  Iced 


!-eedv  leva.  nr  t of  a protot  *pe  evlce  ro  move  over  the 


Reflect  toe  hula  wit haw t tt|alfir«tlr  laumftlat  traffic 
mil  irutly  faellitat*  the  la  format  law  gat  ha  ring  phaaa 
TM»  mil  ilailUMouh  imtll*  • c«hn  li— 1 mater  or 
itwlarl  wMcl*. 

4.  teuirrh  dmull  ha  llrtctal  towards  dataralalag  tha 

lafloaacaa  of  tha  coapaaaata  of  a pawaaaat  awkgrada  grata* 

oa  tha  k , c,  sad  a oaramtsrs  Tha  aaccaaa  of  this  aadaaror 
wawld  ba  a it#*  towards  opt  1*1 1 lag  daaiga  procadwraa  aad 
tha  jm  ef  paring  aatarlals. 
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•CASUKMCVT  ST  ST  DC 

r»  ftrtalbMrl  noda  1 of  the  LCD  iritn  was  euhcoat  r or  tod  by 
Nr4w  Uilmilty  to  Setose*  Applications  Incorporated  (SAI). 

<ktleo,  during  the  period  'toy,  IS?*  to  OteMtor,  197V  The  root  oats 
of  Section  • of  this  Apytodli  1*  based  principally  ©a  SAI*  final 

ftapert  r I* f arose*  H* 

During  t ha  period  !a**an  1973  to  lose  1974,  the  breadboard  node 
oloct rente*  upgr odod  at  furdue  Ouiversity  and  war*  wa*  done  to 

convert  the  static  naoaera— nt  system  la  order  to  naka  It  nor* 
ce^et  1 b 1 • with  a nay  lap  Sean- no  ring  load  concept  Soctioa  b of  this 
Append  1 1 prooente  a *«n*n  of  this  work. 

Section  c of  thta  Appendix  present  * none  details  of  the  Linear 
Variable  Differential  Transformer  Systen  (LVDT)  developed  at  Furdue 
Uniwerslty  by  G.  baladl  deference  32)  for  a study  sponsored  by  the 

• The  UP  Systen 

Several  optical  systens  war*  considered  during  t hr  development 
ef  the  LCD  systen  These  Include 

1 A retrof lector  nodule  whose  mttlon  1*  detected  by  s laser 
Sen 

2 A United  version  of  the  global  Moire  method 

3 focal  point  shift  detection 

* Trlengulat ion  o'  s light  spot  on  s surface 
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As  tm  trt«pilati«i  atM4  «m  finally  m1mu4  f*t  utilw— t , 
only  tM  iMHn  4*t«lli  co*c<raii|  till*  MtM  ara  pmait*4  Mttla 
Tha  —4a  1 Mfletlai  tha  trl— cnlatlaa  — tho4  la  ib«n 

In  Figaro  A-l  This  ar  rang  a— at  has  — laagi  shift  d^  (Figaro  A-2) 
ptrvaadlnlar  to  tha  11—  of  via*  which  la  ralata4  to  tha  pavo— at 
4a f o r— t 1 on  fey  tea  aquation, 

2 * ala  > tan  <M  - y) 

wharo  a la  tha  — galf  1 cat  lee  of  tha  4a  format  loo  Hhar  t— on  or** 
designated  In  Figaro  A-2.  The  pavaaaat  roughness  fanctloe  la  repre- 
sented hv  f(x) 

Tha  straight  11—  of  tha  ligat  caatar  ray  la  given  by 
y - (»-h)  tan  (40  - y) 


y • a tan  (40  - y)  ♦ a 

THa  lntcraactlon  of  tha  Input  rat  with  f (*)  la  at  giving 

f(*  ) - a tan  (40  y)  ♦ a 

aw!  tha  lntaraeetlon  with  tha  llaplacad  f(a)  la  at  * , giving 

i ♦ ffa^)  * tars  (40  - y)  ♦ a 

alnr  a Taylor  expansion  naar  a to  gat  f(a.)  and  dropping  hlghar 
orders  glvao 
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Hx  ) -'<*)*  — r <X  X ) ♦ <X  ~ X )* 

I O l AO  * i o 

ubstltute  this  Into  the  »e cone  inter sect too  equation  «k  i^tn  t inj 
t »«  tw*  jiwi 

r«*  > 3 

(x  ~ X ) tan  (*>  - v)  - d • f <x  )(x  x ) • — (x,  - x } 
lo  • 1 o 2 ~ e 


which  Is  solved  to  ofctaln 


(x  x ) - 

1 e 
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f (*  ) t*w'*0~>)  <■  f (x  tw'.!^  x>  ' - 2f  («  )l 

o—o  e 


The  tune  point  Is  eh© f ted  fee  the  • opponent  , iwrw,  to  the  axis  of 
t he  In.  vine  lens  where. 


i . * (4I  „ 

2 sin  > *1  o 


"he  i*4|t*  ao-  ton  car,  then  he  obtained  a«  a fuc'ctl on  ot  the 

arraapenen:  wrtnrctr  » ' , the  surface  f(x  and  the  surface  displace 

© 

sent  (}  he  the  equation 

4 - a 4 *-tn 

2 i sin  > 

Haweeer  the  aeeenent  is  reuch  on  the  Sit  of  the  ‘ocu.sed  spot  oe  th* 
rutswar  Tc  •lin.lnate  the  depeodeo  * or  surface  row* rasas  there  are 
tw<  alternative*  First  n»«e  the  fncueed  spot  a focuee  lne  leap 

enentti  that  the  effects  of  rcw'xeu  are  aserared  t xero  and  secoaf 
reduce  the  artp*e  between  the  paeenent  oem*  mtv*  Input  ra*  tc  *ero 
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The  sensitivity  of  the  image  motion  is  maximized  when  3 goes  to  90°  but 
that  is  not  physically  possible.  Also  the  collection  of  light  for 
imaging  depends  on  the  solid  angle  of  the  imaging  lens  with  respect  to 
the  pavement  spot.  The  larger  3 forces  the  solid  angle  to  smaller 

values.  i 

A triangulation  arrangement  using  a laser  input  was  breadboarded 
and  tested.  Figure  A-3  shows  the  laser  setup  in  the  laboratory  at  SAI. 

Some  typical  displacement  signal  curves  are  given  in  Figure  A-4.  These 
data  are  for  four  different  samples  and  y • 23°,  3 * 45°.  Several 
other  arrangements  and  the  theory  given  above  demonstrated  that  these 
calibration  curves  become  smooth  and  repeatable  when  y ■ 0°,  3 * 45° 
using  a cylinder  lens  to  focus  a line  on  the  surface.  Also  shown  in 
Figure  A-5  is  a curve  which  was  obtained  by  using  the  sum  of  the  two 
position  signals  as  a normalization  factor.  An  extension  of  the  linear 
range  is  obtained  with  the  normalization.  Using  the  normalization 
method  also  yields  a further  more  Important  advantage,  that  is,  that 
the  reflectivity  of  the  surface  is  removed  from  affecting  the  signal. 

Figure  A-6  shows  normalized  data  for  concrete  and  gray  sandpaper  which 
have  the  maximum  and  minimum  reflectance  of  the  sample  set  used. 

Their  reflectance  difference  was  greater  than  a factor  of  2.  The 
differences  in  the  two  signals  appear  in  the  regions  where  normalization 
has  extended  the  linear  range.  The  major  contributor  to  this  error 
was  the  Inaccurate  measurement  of  the  normalization  constant.  If  higher 
precision  electronics  had  been  used  in  this  breadboard,  most  of  the 
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Figure  A— 5.  Normal  Incidence  Displacement  Curve 


signal  differences  would  disappear 
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The  choice  was  made  to  use  LED  sources  in  a triangulation 
arrangement  to  deflections.  This  arrangement  is  shown  in  Figure  A- 7 
and  is  the  same  triangulation  arrangement  as  was  tested  before  except 
that  the  light  emitting  diode  replaces  the  input  laser  beam  and  the 
cylinder  lens  has  been  modified  by  adding  a convex-piano  spherical 
lens  to  form  a composite  spherical-cylinder  lens.  This  lens 
arrangement  was  developed  by  mapping  the  three  dimensional  distribution 
of  energy  using  an  aperture  and  photodiode.  The  peak  of  the 
intensity  was  found  for  each  lens  arrangement  until  an  optimum 
arrangement  was  determined.  The  displacement  signal,  shown  in 
Figure  A-8  demonstrates  effects  of  the  change  in  sharpness  of  the 
focal  zone.  The  curve  shows  a change  in  slope  in  either  side  of  the 
central  focus  region  that  is  due  to  the  larger  area  for  the  "line" 
focussed  on  the  surface. 

The  basic  optical  principle  as  described  herein  remained 
unchanged  throughout  the  course  of  system  development.  The  electronics 
and  the  mechanical  attachments  were  however  modified  to  be  more 
consistent  with  the  concepts  of  the  moving  load-moving  beam  which 
evolved  subsequent  to  the  conclusions  of  SAl's  effort, 
b.  Additional  Work  on  LED  System 

Repackaging  of  the  breadboard  electronics  and  modification  of  the 
software  supplied  by  SAI  was  done  by  personnel  at  Purdue  University. 
This  work  consisted  of  the  following: 

1.  Converting  the  basic  wiring  of  the  breadboard  electronics 
to  printed  circuits.  Figure  A-9  shows  the  repackaged 
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Figure  A— 7.  Light  Emitting  Diode  Triangulation  Arrangement 


Figure  A-8.  Light  Emitting  Diode  Displacement  Curve 
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electronics. 

2.  Testing  the  basic  moving  beam-moving  load  concept. 

Figure  A-10  shows  the  basic  test  setup. 

3.  Changing  the  software  to  be  compatible  with  the  new 
concept  of  a moving  beam-moving  load.  Figure  A-ll  shows  a 
flow  chart  of  the  modified  software. 

c.  The  LVDT  System 

The  LVDT  system  was  designed  and  fabricated  at  Purdue  University 
as  part  of  a study  sponsored  by  the  Federal  Highway  Administration 
(FHWA).  Because  an  overview  of  the  system  has  already  been  presented 
previously  in  this  repot  only  the  details  of  the  system  components 
are  present  here. 

The  use  of  the  system  to  measure  response  functions  has  been 
discussed  previously  in  this  report.  For  additional  details  see 
Baladi  (Reference  52) . A listing  of  the  basic  components  of  the  LVDT 
system  is  as  follows: 

•Light  Weight  Aluminum  Z-beam,  80  inches  long,  12  inches  wide 
•Linear  Variable  Differential  Transformers  (LVDTs) 

•2  Way  Screw  Jack 
•Pipe  (pivot) 

•Small  Screw  Jacks  under  LVDTs 
•6-Channel  Oscillograph  Light  Beam  Recorder 
•Kodak-2022-light  Sensitive  Paper 
•Sears  Electric  Power  Generator 
•AC- DC  Power  Converter 


• Micrometer-Caliper 


Figure  A-10. 


Test  Setup  Moving  Beam  - Moving  Load  Operation: 
Purdue  University 
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Figure  A-ll.  Summary  Flow  Chart  of  LED  System  Software 
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APPENDIX  B 

SERIES  A:  FIELD  INVESTIGATIONS 

The  Series  A field  investigation  program  consisted  of  conducting 
field  tests  at  Eglin  Air  Force  Base,  Florida  and  Pease  Air  Force  Base 
in  New  Hampshire  during  the  period  of  this  study. 

a.  Eglin  Air  Force  Base 

Two  sites  were  selected  at  Eglin  Air  Force  Base.  The  location  of 
the  sites  is  shown  in  Figures  B-l(a)  and  B-3(a).  At  each  site,  LVDT 
gages  were  embedded  in  the  pavement.  At  Site  1,  five  gages  were 
installed  and  at  site  2,  four  gages  were  placed  in  the  pavement.  The 
locations  of  the  gages  are  shown  in  Figures  B-l(b)  and  B-3(b). 

Details  of  the  boring  locations  and  site  profile  characteristics 
are  shown  in  Figures  B-2  and  B-4. 

A summary  of  the  gage  installation  procedure  is  shown  in  the 
photosequence  in  Figure  B-5.  Figure  B-6  shows  the  details  of  typical 
Installed  gage.  After  Installation,  the  gages  were  calibrated  as 
discussed  in  Boyer  (Reference  11) , and  pavement  response  functions 
v ire  obtained  as  discussed  previously  in  this  report  (see  section  5.1). 

Grids  were  laid  out  as  shown  in  Figures  B-l(b)  and  B-3(b)  and 
response  functions  were  determined  using  the  LVDT  beam.  The  gages  were 
first  calibrated  using  the  equipment  shown  in  Figure  B-7 . Then 
response  functions  were  obtained  at  each  of  the  grid  locations  shown 
in  Figures  B*»l(b)  and  B-3(b).  Figure  B-8  shows  a typical  pass 
being  made  by  the  P-2  fire  truck. 
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Figure  B-l.  Eglin  Air  Force  Base:  Site 


Figure  B-3.'  Eglin  Air  Force  Base:  Site 


Details  of  Boring  Location  and  Site  Profile  Characteristics 


(a) Coring  Pavement 


(b) Extract  Cores 


Figure  B-5.  Gage  Installation  Procedure 


(<J)  Paraphernalia  Ready  for  Installation  in  Hole 


Figure  B-5.  (Continued) 


(g)  Backfill 
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(j)  Backfill  Around  Upper  Paraphernalia 

Figure  B-5.  (Continued) 
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(n)  Typical  Installation- Partially  Epatcied  Cables 


Figure  B-5.  (Continued) 


(o)  Epoxy  Coble  in  Groove  to  Edge  of  Povement 


(p)  Gages  Ready  for  Element  (Sensor)  Installation 


Figure  B-5.  (Concluded) 


Figure  B-8.  Photo  of  P-2  Fire  Truck  Making  a Pass  During  LVDT 
Beam  Testing 
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At  the  locations  shown  In  Figure  B-l(b)  and  B-3(b),  attempts 
were  also  made  to  estimate  the  pavement  transfer  function  by  dropping 
a 60  pound  weight  onto  the  pavement  close  to  the  gage  at  the  tip  of 
the  beam. 

b.  Pease  Air  Force  Base 

One  site,  on  the  main  runway,  was  tested  at  Pease  AFB  in  New 
Hampshire.  At  this  site  only  the  LVDT  beam  was  used  to  obtain 
measurements  of  pavement  response  functions  in  the  manner  previously 
described,  using  the  P-2  fire  truck  and  the  KC-135A  aircraft  as  prime 
movers.  Typical  deflection  measurements  are  presented  in  Figures  B-9 
and  B-10. 


APPENDIX  C 


COMPUTER  PROGRAMS 

This  Appendix  contains  the  details  of  the  computer  program  for, 

1.  Signature  determination 

2.  Parameter  determination 

3.  Prediction 

a.  Signature  Determination 

The  signature  as  defined  in  the  text  is  calculated  using  the 
equations  (3)  through  (10)  and  the  subroutines  EXACFIT  and  SIGNTR. 
The  procedure  is  as  follows: 

With  reference  to  Figures  C-l  and  C-2,  a value  of 
r eak(F^P)^  *8  aS8ume^  (0*2)  then  the  left  hand  side 
(FLHS)  and  the  right  hand  side  (FRHS)  of  equation  5a 
are  calculated.  Their  difference  TVAL(l)  is  then 
calculated,  is  then  incremented  by  0.0125 

(DFXP)  and  the  left  hand  side  and  right  hand  side  of 
equation  5a,  and  their  difference  TVAL(2)  is  again 
calculated.  The  ratio  (FDIF)  of  TVAL(2)  and  TVAL(l) 
is  calculated.  If  this  ratio  is  positive, 
is  incremented  by  0.0125  (DFXP)  and  the  ratio  FDIF  is 
again  found.  This  procedure  is  repeated  until  the  ratio 
is  first  negative.  This  is  shown  as  S in  Figure  C-2a; 


^"  Symbols  in  capital  letters  correspond  to  symbols  used  in  the  sub- 
routine EXACFIT  which  performs  the  calculations. 
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TVAL(j) 


peak 
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Figure  C-2.  Details  of  Interpolation  in  Subgrade  EXACF1T 
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the  point  immediately  before,  at  which  the  ratio  Is  positive 
Is  labelled  L.  In  Figure  C-2b,  the  point  before  the  ratio 
becomes  negative  is  L'  and  S'  is  the  point  at  which  the 
ratio  is  negative.  The  point  where  the  difference  between 
the  left  hand  side  and  the  right  hand  side  of  equation  5a 
is  zero  is  found  by  linearly  interpolating  between  L*  and 
S'.  With  reference  to  Figure  C-*2b 


TFDIF  - | TVAL  (J-l)  | + | TVAL  (J)  | 


FFXP  , , - FXP  - ‘ DFXP  - r*  . 

required  j TFDIF  peak 


^peak  *s  t^en  calculated  from  equation  5 where 

in  z - in  z 

max.  max. 

a i — ± £. 

peak  r'  . r'  . 

peak  peak 


(5a) 


and,  Is  calculated  from  equation  3a 


A'  , - z /[exp  (8*  . x 


rpeak 


peak  max. 


peak  1 


)] 


(3a) 


This  procedure  is  then  repeated  using  (x^,  z^) , (x^,  z^)  and 
(x5,  z5)  instead  of  (Xj^,  z^) , (x^t  z and  (x^,  z^)  • and  values 


of  r"  , , 3"  , and  A . are  calculated, 

peak  peak  peak 


The  desired  r , , 3 , and  A , are  found  as  the 

peak’  peak  peak 


average  of  the  values  using  the  equations, 


r . “ (r*  . r"  . )/2 

peak  peak  peak 


’ (fW  + 6p.ak)/2 
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Vale  * <Vak  + A;«k)/2 


The  A . , r , , and  6 , values  thus  calculated  are 

peak*  peak'  peak 

the  parameters  describing  the  maximum  lateral  deflection  profile. 


The  signature  is  then  calculated  using  the  r ^ 


parameter  and  the  basic  equation,  equation  2a  but  in  the  form 


A(t)  - z(t)/ exp  (6peak  Xj^  peak) 


(10) 


where  z(t)  refers  to  the  deflection  at  gate  1 (closest  to  the  wheel), 
and  t is  the  time.  Subroutine  SIGNTR  performs  these  computations  by 
direct  substitution  of  z(t),  x^  rpeak  and  ®peak  in  e9umtl°n  (10). 


b.  Parameter  Determination 


The  k,  c and  m parameters  as  defined  by  equation  (1)  in  the 
text  are  calculated  using  the  methodology,  concepts  and  equations 
presented  in  Section  III  in  the  paragraph  with  subheading  "Parameter 
Determination".  The  subroutines  EXACT1,  NWCONV,  SIGNTR,  and  EXACFIT 
are  called  from  the  main  program  FNDPRM  in  order  to  perform  the 
necessary  computations. 

With  reference  to  Figure  C-3  and  Figure  C-4,  initially,  a value 
of  the  spring  constant  k is  assumed  where 


k - EQKA  - AA1  x PMAX/DMAX 


(29) 


,18 


where  AA1  - 0.05 

PMAX  - 25000,  and 

DMAX  * Maximum  value  of  the  signature. 


18 


For  convenience,  the  symbols  in  capitol  letters  in  section  (4)  are 
the  same  symbols  used  in  the  computer  program  FNDPRM. 
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Then  a value  of  m Is  calculated  using  the  equation 

m - EQMA  - AA3  x PMAX/DRV2  (IMAXD)  (30) 

where  DRV2  (IMAXD)  is  the  acceleration  at  maximum  deflection  and  AA3  - 
.91.  Equation  30  is  derived  as  follows:  Let 

k - AA1  x PMAX/DMAX  (29) 

at  maximum  deflection,  i.e. , section  (3),  Figure  7.  Using  a constant 
AA2 , the  equation 

AA2  x PMAX  - EQMAXDRV2 (IMAXD)  + EQKA  x DMAX 
- EQMAXDRV2 (IMAXD)  + AA1  x PMAX 

is  obtained.  Therefore 


AA2  x PMAX  - AA1  x PMAX  - EQMAXDRV2 (IMAXD) 

(AA2  - AA1)PMAX  - EQMA  x DRV2 (IMAXD)  (31) 


Letting,  AA3  ■ AA2  - AA1  yields 


AA3  x PMAX 
DRV  (IMAXD) 


EQMA 


Next  a value  of  c is  calculated  using  equation  (32)  where 


c - EQCA  - EQMA  x DRV2(II)/DRV1(II)  (32) 


where  DRVl(II)  is  the  velocity  at  time  increment  II  (Figure  C-3).  The 
initial  value  of  II  corresponds  to  the  point  of  zero  acceleration  + 1 
(Figure  C-3).  Using  the  initial  value  of  k,  m and  c calculated  using 
equation  (29),  (30)  and  (32)  respectively,  the  transfer  function  is 
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calculated  using  equation  (11a)  or  (lib)  whichever  Is  applicable. 
Subroutine  EXACT1  performs  this  computation. 

Next  the  value  of  (equation  33)  is  calculated  using  the 

discrete  implicit  form  of  the  convolution  integral  equation  (34) 


0(t)  ■ I(t-x)'G(x)  dx 


where 


I(t)  - 


Its  implicit  form  is 


k-i  -1 

I(i)  - 0(i)MT  - l G(m)  *1  (k  + 2 - m) 
nr»l 


where  the  initial  value  of  i is  2,  and  1(1) 


I(i)-AT 


(34a) 


G(x)  is  the  transfer  function  and  AT  is  the  time  increment.  Sub- 
routine NWCONV  performs  these  calculations. 

Next  F(t),  the  equivalent  input  function,  is  calculated  using 
the  equation 


F(t)  - m x I(t)  - DF(I) 


(34a) 


The  maximum  value  CPMAX  of  F(t),  is  then  found  and  compared  with 
25,000  lbs.  If 


ABERFG  - CPMAX  - 25000  > 500 


then  II  is  incremented  by  1 , and  the  procedure  subsequent  to  that  step 
is  repeated  using  the  DRV2(II)  and  DRVl(II)  values  corresponding  to 


fe&OK&v-: 


the  new  II.  ILP,  JJ  and  II  are  counters  which  have  maximum  values  of 
40,  20  and  IMAXC  (Figure  C-3)  respectively.  These  counters  ensure 
that  the  tested  range  of  values  for  AA1  is  from  .05  to  1.05,  for  AA3 
it  is  from  -.91  to  -.01,  and  for  II  it  is  from  the  first  point  after 
the  point  of  inflection  to  the  maximum  deflection  (Figure  C-3).  If 

ABERFG  - |CPMAX-25000|  < 500  (35b) 

then  c is  calculated  using  the  equation  for  section  (1),  Figure  9, 
equation  36 


CMAXV  - (ZVALUF  - EQKA*  ZVALV0)/ZVALU1  (36) 


where  CMAXV  *•  c 

ZVALUF  * force  at  section  (1)  (Figure  7) 

EQKA  - k 

ZVALVO  ■ the  deflection  at  section  (1)  (Figure  7) 

ZVALU1  ■ the  velocity  at  section  (1)  (Figure  7) 

The  ratio  given  by  equation  (37)  is  next  calculated 

CDIFFR  - (CMAXV  - EQCA) /CMAXV  (37) 


where  EQCA  is  the  c value  computed  in  equation  (32).  For  the  first 
value  of  CDIFFR,  the  procedure  increments  II  by  1 and  repeats  all 
the  steps  from  that  point.  For  the  second  and  subsequent  values  of 
CDIFFR,  the  ratio  of  the  current  and  the  preceding  value  is 
calculated  using  the  equation 


CRATI0  - QC(ICC)/QC(ICC-1)  (38) 
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where  QC(ICC)  refers  to  the  current  ratio  of  CDIFFR,  AND  QC(ICC  - 1) 
refers  to  the  preceding  ratio  of  CDIFFR.  If 

CRATIO  > 0 

then  the  difference  between  the  c values  computed  from  two  Independent 
equations,  equation  (32)  at  the  crossing  point,  and  equation  (36)  at 
the  point  of  Inflection,  has  not  passed  through  zero,  so  II  Is  Incre- 
mented and  the  subsequent  steps  are  repeated.  If 

CRATIO  < 0 

then  the  difference  between  the  c values  calculated  by  equations  (32) 
and  (36)  has  passed  through  zero  and  the  c values  calculated  using 
these  two  equations  are  approximately  the  same. 

It  has  been  found  that  the  two  controls,  namely, 

ABERFG  < 500,  and, 

CRATIO  < 0 

ensures  that  the  conditions  stated  In  paragraph  3.2.2,  are  satisfied. 

Finally,  the  exact  solution  for  the  force  is  calculated  by 
substitution  EQCA,  EQKA,  EQMA,  DRVO,  DRV1,  and  DRV2  in  the  equation, 

CMPXX(I)  - EQKA  x DRVO(I)  x EQMC  x DRV1(I)  + EQMA  x DRV2(I)  (39) 

where  CPMXX(I)  is  the  force  at  each  time  Increment  I In  the  record. 

The  force  at  the  crossing  point  calculated  by  the  convolution  and 
substitution  technique  are  computed.  A listing  of  the  program  FNDPRM 
and  a typical  output  Is  given  at  the  end  of  this  Appendix. 

220 


sarARfJuai 


nvmir.:.  j.  wj lvw w ettwi i?  ■'>  ; ..•  ..  v_v.v. 


c.  Prediction 

The  subroutines  EXACFT,  SIGNTR,  NWCONV,  EXACT  1,  ERECT2,  and 
MODPRM  are  used  by  the  main  program  PREDCT  to  predict  the  signature 
of  a vehicle  (Vehicle  A)  at  a site  (site  I)  over  which  the  vehicle 
has  not  traversed.  It  is  required  to  know  the  signatures  of  vehicle 
A at  another  site  (standard  site),  and  a standard  vehicle  (vehicle  B) 
at  both  the  standard  site  and  site  I.  Since  details  of  the  methodology 
for  determining  the  parameters  of  the  transfer  function,  the  loading 
function,  and  the  equivalency  function  has  been  detailed  previously 
only  a summary  flow  chart  of  program  PREDCT  is  given  here.  Figure 
C-5  shows  the  summary  flow  chart  of  program  PREDCT.  A listing  of  the 
program  and  a typical  output  is  given  in  the  following  pages  of  this 
Appendix. 
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Figure  05.  Sumery  Flow  Chert  of  Program  PREDCT 
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PROGRAM  FNDPRM< INPUT. DUTPUT. TAPES* INPUT. TAPE6=0UTPUT» PLOT) 


THIS  PROGRAM  COMPUTES  PARAMETERS  OF  THE 
REDUCED  TRANSFER  FUNCTION 
FROM  DEFLECTIONS  MEASURED  AT  5 POSITIONS 
PERPENDICULAR  TO  THE  DIRECTION  OF  TRAVEL 
OF  THE  LOAD  VEHICLE 


SUBROUTINES  REQUIRED 

<1>  EXACT1 
<2>  NWCONV 
<3>  EXACFIT 
<4)  SIGNTR 


INPUT- 


NR— NUMBER  OF  RECORDS 
JM — NUMBER  OF  GAGES 

AIN— DISTANCE  OF  GAGE  FROM  OUTER  EDGE  OF  TIRE 

NP-NUMBER  OF  DATA  POINTS 

PT— TIME  INCREMENT  IN  SECONDS 

UMAX — MAXIMUM  WHEEL  LOAD  USED 

INITX— DIGITIZED  DATA 

CALG— CALIBRATION  DATA 


PR 


2. 


3. 

4. 

5. 

6. 

7. 

8. 
9. 

10. 

11. 

12. 

13. 

14. 


DIMENSION  FRESP<11» 100.5). INITX<1 1.1 00. 5>.CALG< 11.5). NP< 11). PT<11> 
1 . UMAX< 1 1 ) »DRV0< 1 00) > GMAX< 11.5) 

DIMENSION  VIPY<200> »VISIG<200) 

DIMENSION  AI <5>»  FI 
DIMENSION  AIN<5.11> 

DIMENSION  QC<2000> 

DIMENSION  APKClOOJ.BTClOOXTXPONdOO) 

DIMENSION  FR< 1 0) . DFY<  1 0) . FRR< 10.6) 

DIMENSION  DF<200> 

DIMENSION  T <200) 

DIMENSION  0<200>,GG<200> 

DIMENSION  V<200>.G<200> 

DIMENSION  IVEL<11»5> 


13. 

16. 

17. 

18. 

19. 

20. 
21. 
22. 

23. 

24. 
23. 
26. 

27. 

28. 
29. 


C 

C 


DIMENSION  DR VI <200 J .DRVE^POO). CPMXX<50) 


READ<3. 804)NR»  JM 
DO  810  K=1 » NR 

READ<5.29XAIN<J.K>.  J=l.  JM) 
READ<5. 801 >NPOO . PT<K> » WMAX<K) 
JP=NP<K) 

DO  800  1=1. JP 

READ<5.802XINITX<K.  I. J>. J=1.S> 

800  CONTINUE 

READ<3.803XCALG<K.  J>.J=1.5> 

810  CONTINUE 

801  FORMAT  < 15. 2F1 0. 4> 

802  FORMAT <b<I4.2X>> 

803  FORMAT <5F12.10> 

804  FORMAT<2I10> 

29  FORMAT <6F10.5> 
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C####*##*####*RECORDS  RUN 
C 

C EGLIN  SITE  1 LOC  B-4  LD.  CT.  3.10.76  3.0  .2 

C 

Q3K)K9K3li)K)KiKil3K)li3IE3KiKSK 

Q3K)K)K)KiK)KIK)K3t()Kil3i()K)K 

c 

TOP  OF  MAIN  LOOP  THAT  OPERATES  ON  ONE  RECORD  AT  A TIME 


C 

C 


C 

C 

C 

C 

C 


DO  811  K=1.NR 


CALCULATES  DEFLECTION  RESPONSE  FUNCTION  FOR  JM  GAGES  AND 


THE  PEAK  FOR  EACH  GAGE  FOR  THE  K TH  RECORD 


nn  RflA  .1=1  ..IM 
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■u\  i 


r\ 

% 


m 


,.n.J)=AIN<J,K> 

RMAX=0.0 

JP=NPO<> 

DC1  807  1=2, JP 

FRESPCK,  I-l»  J)=FLOAT<INlTX<K,I»  J)-INITX<K» 1,  J))*CALG<K»  JV12. 
VIPY<I-1)=FRESP<K,  1-1,2) 

IF<FRESP<K» 1-1, J).LT.RMAX)GO  TO  807 
RMAX=FRESP<K,I-1,J> 

807  CONTINUE  * 

GMAXCK, J)=RMAX 
FI < J)=GMAX<K, J) 

806  CONTINUE 

LISTS  REDUCED  DATA 

WRITE<6»  833) 

833  FORMAT < 1H1 » 40X, *NEW  RECORD*/'/') 

URITE<6,  814)' 

814  FORMAT </"^lX»*NUMBER  OF  DATA  POINTS  RECORD  NUMBER  MAX.  UH 

1EEL  LOAD  TIME  INCREMENT  NUMBER  OF  GAGES*/') 

URITE<6,815)NP<K>,K,WMAX<K),PT<K), JM 

815  FORMAT < 1 OX, 14, 19X, 14. 18X, F6. 0, 19X> F4. 2, 18X, 12/-) 

URITE<6,  816) 

816  FORMAT </^4 OX, *DEFLECT IONS  IN  FT*/-) 

URITE<6»  817) 

817  FORMAT < 13X,  *GAGE  1 GAGE  2 GAGE  3 GAGE  4 

1 GAGE  5*/) 

DO  818  1=1. JP 

WRITE<6, 819) <FRESP<Kj I, J), J=1 » JM) 

818  CONTINUE 

819  FORMAT<10X,5<F7.5, 10X>> 

URITE<6»  820) 

820  FORMAT </'/'4 OX, *PEAK  DEFLECTIONS*/') 

URITE<6»819>  <GMAX<K, J), J=1 , JM) 

FIT  MEASURED  PEAK  LATERAL  DEFLECTION  PROFILE  TO 
EQUATION  USING  SUBROUTINE  EXACFIT 


CALL  EXACF I T<F I , TFXP, TBETAF , FAF , AI ) 

APK<K>=FAF 
BT<K)=TBETAF 
TXPUN  < K ) =TFXP 

VBETA=TBETAF 

VRVAL=TFXP 

VIPX=AIN<£»K) 

NV=JP 

CALCULATE  SIGNATURE  OF  PROTOTYPE  USING 
SUBROUTINE  SIGNTR 


CALL  S I GNTR < VBETA , VRVAL .VIPX.VIPY.VISIG.NV) 


CALCULATE  FIRST  AND  SECOND  DERIVATIVES 


DO  74  1=1, JP 
DRVO<I)=VISIG<I) 
TRMl=<VISIG<I)+VISIG<I  + l))/'2. 
IF<I.EQ.  1 )1RM£=VISIG<  I )/'2. 
TRM2=<VISIG<I)+VI.SIG<I-l))/'2. 
DRV1  < I )=<TRM1-TRM2)/'PT<K) 

74  CONTINUE 

DO  75  1=1, JP 

TRN1=<  DRV1 < 1+1 )+DRVl < I))/2. 
IFCI.EQ.  l)TRN2=DRVl<I>/'2. 
TR.N2=  < DRV 1 < I ) +DRV  1<  1-1)  )/£. 

DR  V2  < I ) = < TRN  l-TRNO/'PT  <K) 

75  CONTINUE 


IDENTIFY  FIRST  PEAK  BY  FINDING  POINT  WHERE  FIRST 
DERIVATIVE  CHANGES  SIGN 
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c 


1 


BOTM=0. 00000000001 
DO  726  1=1, JP 
HUMP1=DRV0<I+1>-DRV0< I ) 
IF<HUMP1.GE. 0.  0>GO  TO  726 
JPP=I+6 
GO  TO  727 

726  CONTINUE 

727  CONTINUE 


REPLACE  ZEROES  BY  AVERAGING  BEFORE  AND  AFTER  POINTS 


91. 

92. 

93. 

94. 
CAUTION 

95. 
CAUTION 

96. 

97. 

98. 

99. 
100. 
101. 


I MAXV=  0 
DV1MX=0. 0 
DO  315  1=1, JPP 

IF<DRV1<I).EQ.0.  0>DRVKI>=<DRVl<I  + l>+DRVl<I-l>>/'2. 

IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  ******************************* 

IF <DRV1 <1 > . EQ.  0.  0)DRV1 <I)=<DRV1 < I+2>+DRVl  < I-2JV2. 

IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  KnKimminiiiKiininwiinwiinwwiKiKiKituwauitwMwiinwmmiiniiw 

IF<ABS<DRV1<  I > >.LT. BOTH) DRV 1< I >=<DRVKI+l>+DRVl<I-l>>/'2. 

IF<ABS<DRV2<I > > . LT. B0TM)DRV2< I >=<DRV2< 1+1 >+DRV2< 1-1 > V2. 

IF<DRV1<I>.LT.DV1MX>G0  TO  315 

DV1MX=DRV1<I> 

IMAXV=I 
315  CONTINUE 


FIND  MAXIMUM  VALUE  OF  SIGNATURE 


DMAX=DRV0< 1 ) 

DO  53  1=2* JPP 

IF<DRVO<I>.LT.DMAX)GO  TO  53 
DMAX=DRVO<I > 

IMAXD=I 
53  CONTINUE 


REPLACE  ERRONEOUS  ACCELERATIONS  BY  AVERAGING 
BEFORE  AND  AFTER  POINTS 


108. 

109. 

no. 

111. 

112. 

113. 

114. 

115. 

116. 

117. 

118. 

119. 
CAUTION 

120. 
CAUTION 

121. 


I MAXV 1=1 MAX V- 1 
DO  57  1=2, IMAXV1 

IF<DRV2< I > . LE. 0. 0)DRV2< I )=<DRV2<I-1 >+DRV2< 1+1 >>/2. 

IF<IiRV2< I ) . LE , 0. 0>DRV2< I >=<  DRV2< I -2)+DRV2< I+2>  >/2. 

57  CONTINUE 

I M AXVA= I MAXV+ 1 
DO  109  I=IMAXVA, JPP 

IF<DRV2< I ) . GE. 0. 0)DRV2< I ) = <DRV2< 1-1 >+DRV2< 1 + 1 > V2. 

I F < DRV2 <I).GE.  0.  0) DRV2 < I ) ■ < DRV2 < I -2 > +DRV2 <I+2>V2. 

109  CONTINUE 

DO  107  1=1, JPP 

IF<DRV1 < 1+1 ) . EQ. DRV1 < I ) )DRV1 < 1+1 )=<DRV1 <I+2>+DRVl < !>)/£. 

IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  mmMmmmmmmmmmmmmmmmmmMmMMmmmmmmm 

IF  <DRV2< I + l ) . EQ. DRV2< I ) )DRV2< 1+1 >=<DRV2< I+2)+DRV2< I ) )/2. 

IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  *m***m*mm*mmmmm*mmmmmmmmmmmmmmm 

107  CONTINUE 


FIND  TIME  INTERVAL  WHERE  VELOCITY  CHANGES  SIGN 
AND  FRACTION  OF  VELOCITYTO  ZERO  BY 
LINEAR  INTERPOLATION 


DO  789  IA=1 » JPP 
DEROT  S=DRV  1 < I A+l  > /”BRV  1 < I A> 
IF<DEROTS.GT. 0. 0>GO  TO  789 
DZERO=DRV 1 < I A > -DRV 1 < I A+ 1 > 
DDTZ=DRV 1 < I A ) /DZERO 
NAFNL=IA 
GO  TO  788 
789  CONTINUE 
738  CONTINUE 


FIND  INTERVAL  AND  TIME  UHERE  ACCELERATION  IS 
ZERO  AND  THE  CORRESPONDING  DISPLACEMENT 
AND  VELOCITY 


DO  389  IA=1» JPP 
ZEROTS=DRV2<  IA+1  >/'DRV2<  IAJ_ 


m. 


\VVj 

r>7.; 

* ; v 

mmm 


‘VUVMWmTHIl1 


tl.  i-.  t..  ir.  I*-  i.  I 


IF<ZEROTS.GT;0.0)GU  TO  889 
TZER0=DRV2 < I 8 > -DRV2< I A+ 1 > 

FDTZ=DRV2<  IAVTZERO 
IAFNL=IA 

DR0RNG=DRV0<IA+1>-DRV0<IA> 
ZVALUO=DRVO<IA>+FDTZ*DRORNG 
DR 1 RNG=DRV 1 < I A+ 1 > -DRV 1 <1 A > 
ZVALU1=DRV1<IA>+FDTZ*DR1RNG 
GO  TO  888 
889  CONTINUE 
888  CONTINUE 
C 

CMKXKKKXXK 

WRITE<6»99> 

99  FORMAT (✓37X.  *FOR  SIGNATURE*^ 

URITE<6»  76> 

76  FORMAT <^7X»  * INCREMENT  SIGNATURE 

lERATION*/') 

MRITE<6» 77>  < I . DRV0< I > . DRV1C I > . DRV2< I > > 1*1 . JPP> 

77  F0RMAT<7X» I14.3E13.5) 

WRITE<6» 1091) 

1 091  FORMAT  < s/s/sssss/ IX » 

IMAXF=0 

PMAX=UMAX<K> 

AA1=0. 05 
ICC=0 


VELOCITY 


ACCEL 


TOP  OF  LOOP  INCREMENTING  AA1 


DO  104  ILP=1»40 
EQKA= A A 1 XPMAX-'DMAX 
AA3—0.91 


TOP  OF  LOOP  INCREMENTING  AA3 


C 

c 

c 

c 


DO  103  JJ=1>  20 
IT*IMAXV 

EQMA=AA3*PMAXxDRV2< IMAXD) 

CPMAX=0. 0 

AA2=AA1+AA3 

IF<AA2.LT. O.OJGO  TO  1021 

IMAXC=IMAXD 

IT*IT+1 


TOP  OF  LOOP  VARYING  TIME  BETWEEN  ZERO  ACCELERATION 
AND  ZERO  VELOCITY  AND  CALCULATING  EQCA 
FOR  EACH  INCREMENT 


DO  102  II*IT. IMAXC 
EQCA*EQMA*DRV2< I I VDRV1 < I I > 
EQCA—EQCA 


EQVK=EQKA 

EQMASS=EQMA 

EQVC=EQCA 

MM=JPP 

DT=PT<K> 


CALCULATING  REDUCED  TRANSFER  FUNCTION  USING 
SUBROUTINE  EXACT 1 


CALL  EXACTKMM.DT.  EQVK. EQMASSi T»  GG.EQVC) 

DO  73  1=1. JPP 
0<I>=DRVO<I> 

V<I)=GG<n 
73  CONTINUE 
N=JPP 

CALCULATE  XT)  USING  SUBROUTINE  NWCONV 
BY  IMPLICIT  CONVOLUTION 


vi  J 


I 


18f.  CALL  NWC0NV<M*7.G»a*I>T> 

c 

c 

C CALCULATING  EQUIVALENT  FORCE AS  PRODUCT  OF 

C I<T)  AND  EQMASS 

c — — — — — 

182.  DO  308  MP=1*JPP 

183.  DF  <MP>=G<MP>*EQMASS 

184.  308  CONTINUE 

185.  DO  509  MA=1» JPP 

186.  DF<MA>»<DF<MA>*DF<MA+l>>/2. 

187.  509  CONTINUE 

C 

C AVERAGING  AND  FINDING  MAXIMUM  EQUIVALENT  FORCE 

C 

188.  CPMAX=0. 

189.  DO  501  M=1 * JPP 

190.  DF<M>*<DF<M>*DF<M+1>>'2. 

191.  IF<DF<M).GT.CPMAX>GO  TO  20 

192.  GO  TO  502 

193.  20  CPMAX=DF<M> 

194.  ICF=M 

195.  502  CONTINUE 

196.  501  CONTINUE 

C . 

C EXECUTING  FIRST  CHECK  TO  ENSURE  THAT  THE 

C ABSOLUTE  DIFFERENCE  BETWEEN  THE  MAXIMUM  COMPUTED 

C FORCE  AND  23  KIPS  IS  LESS  THAN  300  LBS 

C 

197.  ERFCE=CPMAX-PMAX 

198.  ABERFG=ABS  < ERFCE  > 

199.  IF<ABERFG.GT.500. 0>GQ  TO  102 

C 

C EXECUTING  SECOND  CHECK  TO  ENSURE  THAT  THE  DIFFERENCE 

C BETWEEN  THE  C VALUES  COMPUTED  AT  THE  CROSSING  POINT 

C AND  THE  POINT  OF  INFLECTION  GOES  THROUGH  ZERO 

C 

200.  icc*ia>i 

201.  DFRNG*DF  < I AFNL+ 1 > -DF  < I AFNL  > 

202.  ZVALUF^DT  < IAFNL>+FDTZ*DFRNG 

203.  CMAXV- < ZVALUF -EQKA*ZVALU 0 > rZVALU 1 

204.  CD  I FFR*  < CMAXV-EQCA  > -^CMAXV 

205.  ABCD*CDIFFR 

206.  QC< ICC>*CDIFFR 

207.  IF<ICC.EQ.l)GO  TO  102 

208.  CRATIO*QC< ICC)^QC< ICC-1 > 

209.  IF<CRATIO.GT .0. 0)GO  TO  102 

C 

C CALCULATING  FORCE  RATIOS  COMPONENTS  OF  FORCE 

C AND  TOTAL  FORCE  AT  THE  CROSSING.  POINT 

C 

210.  FCHK«EQKA*DRVO<II> 

211.  D I FCHK=FCHK-DF  < 1 1 > 

212.  DFCHKR=DIFCHK/FCHK 

213.  FPFRAC=DKII>/PMAX 

214.  ICDPT*II 

213.  TRMM=EQMA*DRV2  < 1 1 > 

216.  TRMC=EQCA*DRV1 < I I > 

217.  TRMK=EQKA*DRV 0 <1 1 > 

218.  CPMXX< 1 1 >=TRMM+TRMC+TRMK 

219.  WRITE<6»847) 

220.  847  FQR^9T<//'10X,«GOOD  RUN"//) 

221.  GO  Tli  846 
C»****»**»* 

222.  102  CONTINUE 

223.  1021  CONTINUE 

224.  AA3«AA3+ 0.043 

QK**W***K**»«N 

225.  103  CONTINUE 

226.  AAl=AAl+0. 025 

. c******* 

227.  104  CONTINUE 

228.  1000  FQRMAT</'1X,*000*> 

229.  845  F0RMAT<6F15.4» 10X»2I5»F10.3> 

230.  URITE<6» 1002) 
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249. 

250. 

251. 

252. 

253. 

254. 

255. 

256. 

257. 

258. 

259. 

260. 
261. 
262. 

263. 

264. 

265. 

266. 

267. 

268. 

269. 

270. 

271. 

272. 

273. 

274. 

275. 

276. 

277. 

278. 

279. 

280. 

281. 

282. 

283. 

284. 


1002  FORMAT </l OX,  "OUTER  LOOP  COMPLETED*/) 

GO  TO  94 
846  CONTINUE 

URITE<6,1043> 

1 043  FORMAT  < ////////// 1 OX i * I NCREMENT* »2X» *M-COMP . * , 8X, *C-COMP . * » 8X,  *K-C 
10MP. *. 8X» "EXACT*, 8X, "CONVOLUTION*. 8X.*TDT  FUNC*. 8X. "SIGNATURE*/) 

C 

C CALCULATING  COMPONENTS  OF  FORCE  DUE  TO  INERTIAL 

C DAMPING  AND  STIFFNESS 

C 

JPN=JPP-3 
DO  106  1=1, JPN 
TMT=EQMA*DRV2< I > 

TCT=£QCA*DRV1 < I > 

TKT=EQKA*DRVO< I > 

CPMXX< I >=TCT+TKT+TMT 

WRITE <6.1001)1, TMT . TCT » TKT . CPMXX  < I > , DF  < I > » GG< I > , DRV  0< I > 

106  CONTINUE 

C 

C CALCULATING  FORCE  RATIOS  AT  SECTIONS  <1>.<2>.<3> 

C 

AFRNG=DF  < NAFNL+ 1 > -DF  < NAFNL  > 

DVALUF=DF < NAFNL  > +DDTZ*AFRNG 

C"**** 

ZVALUR=ZVALUF/PMAX 

£****** 

DVALZD=DVALUF/,PMAX 

£*"*** 

FDR  T 1 0=AB£RF G/PMAX 

C 

C CALCULATING  TIME  RATIOS  AT  SECTIONSCl). <2>* <3) 

C : 

DENOM=FLOAT  < IMAXD) 

DEN0M1=1. 0 

TCRSPT =FLOAT  < 1 1 )/DENOM 
TZROA-FLOAT  < I AFNL ) /DENOM 
C ******* 

WRITE<6, 1041) 

1041  FORMAT <////////// 1 OX . *EQU I V.  MASS*. 1 0X,*EQUIV.  STIFFNESS*. 9X.*EQU I 
IV.  DAMPING*) 

WRITE<6»2041> 

2041  FORMAT <12X.*SLUGS*. 16X. *LBS/FT*.20X.*LBS-SEC/FT*/> 

WR I TE < 6 . 1 042 ) EQMA . EQKA . EQCA 

1042  FORMAT <1  OX, FI 0.2. 15X.F10.2. 12X.F10. 1) 

URITE<6, 1 044) 

1044  FORMAT <//2 OX,*  FORCE  RAT I OS C FORCE  OVER  25000)*/) 

WRITE<6, 1045) 

1045  FORMAT <25X , *SECT I ON< 1 ) SECTI0N<2)  SECTI0N<3>*/> 

WRITE<6, 1046)2VALUR,FPFRAC»DVALZD 

1046  FORMAT  < 23X , 3F 1 0 . 3 ) 

WRITE<6, 1050) 

1050  FORMAT <///2 OX »*  TIME  RATIOSCTIME  OVER  TIME  TO  MAX.  DEF.)*/> 

DEN0M1=1. 

URITE<6, 1045) 

WR I TE  < 6 , 1 046 ) TZROA , TCRSPT , DENOM 1 
WRITE<6, 1060) 

1060  FORMAT <///20X» *C-PARAMETER  DIFFERENCE  RATIO*/) 

WRITE(6, 106DCDIFFR 

1061  FORMAT (3 OX » F7. 3) 

WRITE<6» 1070) 

1070  FORMAT <///20X,*FDRCE  DIFFERENCE  RATIO  AT  THE  CROSSING  POINT*/) 
URITE<6» 1 061 >DFCHKR 
WRITE<6, 1080) 

1080  FORMAT<///20X» "MAXIMUM  FORCE  DIFFERENCE*/) 

WRITE <6, 1 061 >ABERFG 
1001  FORMAT< 10X, I5.7F15. 5) 

C 

811  CONTINUE 
94  CONTINUE 
C 

STOP 

END 

r 


c- 


I 


1.  SUBROUTINE  EXACFIT<FI»TFXP.TBETAF.FAF»AI> 

C 

C THIS  SUBROUTINE  CALCULATES  THE  PARAMETERS  OF  THE 

C MAXIMUM  LATERAL  DEFLECTION  PROFILE  USING 

C THE  DEVELOPED  INTERPOLATION- ITERATION  PROCEDURE 

C 

C INPUT-  FI-MEASURED  DEFLECTIONS  AT  LOCATIONS  LATERAL 

C TO  THE  TIRE  PRINT 

C AI— LATERAL  DISTANCES  FROM  EDGE  OF  TIRE  PRINT 

C CORRESPONDING  TO  FI 

C OUTPUT  TFXP — R PARAMETER 

C TBETAF— BETA  PARAMETER 

C FAF— A PARAMETER 

C 

2.  DIMENSION  FI<6>.TBETA<3>.TVAL<999>> AI<5>»FA<3>»FFXP<5> 

C$  TRACE  SUBSCRIPTS 

CS  TRACE  ARITHMETIC 

3.  M=1 

4.  IFTX=0 

3.  1=2 

6.  IG=I+i 

7.  IM=I+2 

8.  ITRID=0 
C 


£XXXXXXXXXXXXXXXXXXXXXMXXXXXX 

CXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

c 

9.  710  CONTINUE 

10.  FXP=0.2001 

11.  DFXP=0. 0125 

12.  J=0 

13.  771  CONTINUE 

14.  J=J+1 

15.  FXP=FXP+DFXP 

16.  IFCHK=0 

17.  772  CONTINUE 

18.  776  CONTINUE 

19.  CONSl=ALOG<FI < I > >-ALOG<FI < IG>> 

20.  C0NS2=AL0G<FI < I > >-ALDG<FI < IM> > 

21.  CONF =CONS  1 /'CQNS2 

22.  FLHS=CONF*  < < A I < I )**FXP)-<AI < IM>**FXP>  > 

23.  FRHS=<AI<I)#*FXP>-<AKIG>**FXP> 

24.  TVAL  < J > =FLHS-FRHS 

25.  IF<J.EQ.l)GO  TO  771 

26.  FDIF=TVAL<JVTVAL<J-1> 

27.  ABVAL=ABS<TVAL<J>> 

28.  IF< J. EQ.  999JG0  TO  726 

29.  IF<FDIF.GT. 0. OJGO  TO  771 

30.  TFD I F= ABS  < TVAL  < J- 1 > > + ABS  < TVAL  < J > > 

3 1 . FFXP  < M > =FXP-  < TVAL  < J > /'TFD I F > *DFXP 

32.  773  CONTINUE 
C 

33.  FBETA1=AL0G<FI < I ) )-ALOG<FI < IG>> 

34.  FBETA2=<AKI)X»FFXP<:M>)-<AI<IG?**FFXP<M>> 

35.  TBET  A<M>=FBETA1^FBETA2 

36.  FA1=TBETA<M>*<AI<1>**FFXP<M>> 

37.  FA<M)=FI<  1>/'EXP<FA1> 

38.  GO  TO  777 

39.  IFCITRID.EQ. 1>G0  TO  727 

CAUTION  - STATEMENT  CANNOT  BE  EXECUTED  - NO  STATEMENT  NUMBER  xxxxxxxxxxxxxxxxxxxx 

40.  IF< IFTX.EQ. 1)GD  TO  777 

41.  M=M+1 

42.  IG=IG+1 

43.  IM=IM+2 

44.  IFTX=1 

45.  GO  TO  710 

46.  777  FAF=FA<M> 

47.  TBET  AF =TBET  A < M > 

48.  TFXP=FFXP<M> 

49.  GO  TO  728 

CXXXXKXXXXX 

50.  726  CONTINUE 
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ITRID^I 

M=M+1 

IG=IG+1 

IM=IM+2 

IFTX=1 

IF( IM. LE. 5)G0  TO  730 
URITE<6>  731 ) 

731  FORMAT < /-2  OX » *DAT A DOES  NOT  FIT— CHECK  DATA*/-) 
GO  TO  720 
730  GO  TO  710 

727  CONTINUE 

FAF=FA<2) 

TBETAF=TBETA<2) 

TFXP=FFXP<2) 

728  CONTINUE 
14RITE(6»  775) 

775  FORMAT </-2 OX >*  A PARAMETER  BETA 

UR  I TE  < 6 > 778 )FAF » TBETAF > TFXP 
778  FORMAT (20X»3E15. 5-0 
URITE<6>  750) 

750  FORMAT < ^2 OX > * CALCULATED  Y * MEASURED  Y 
DO  701  M=l*5 
YY 1 =TBE  TAF*  < A I < M ) **TFXP ) 

YY2=EXP<YY1) 

YY=FAF*YY2 

NRITE<6»20)YY»FI<M)> AI <M> 

701  CONTINUE 
720  CONTINUE 
20  FORMAT (20X>  3E15. 5) 


EXPONENT  R*>-) 


LATERAL  DIST.*/) 


RETURN 

END 


SUBROUTINE  SIGNTR<VBETA> VRVAL. VIPX. VIPY. VISIG, NV) 

THIS  SUBROUTINE  CALCULATES  THE  SIGNATURE  BY 
DIRECT  SUBSTITUTION 

INPUT-  VBETA-- BETA  PARAMETER 

VIPX — LATERAL  DISTANCE  OF  GAGE  FROM  OUTER 
EDGE  OF  TIRE  PRINT 

VIPY—NV  VALUES  OF  MEASURED  DEFLECTION 

AT  VIPX  FROM  OUTER  EDGE  OF  TIRE  PRINT 
VRVAL— R PARAMETER 

OUTPUT-  VISIG— NV  VALUES  OF  THE  SIGNATURE 

DIMENSION  VIPY<200)» VISIG<200) 

VTOP= VBE  TA*  < V I PX**VR VAL ) 

DO  40  1=1 »NV 

VISIGO )=VIPY< I >/EXP<VTOP) 

40  CONTINUE 
RETURN 
END 


SUBROUTINE  NUC0NV<N> V.G.O.DT) 

THIS  SUBROUTINE  IMPLICITLY  CONVOLUTES  V<  ) AND  0<  )TO 
YIELD  G<  ) AT  N POINTS 


INPUT-  V— N VALUES  OF  INPUT 
Q— N VALUES  OF  OUTPUT 
DT—TIME  INCREMENT 


OUTPUT- 


G— N VALUES  OF  THE  IMPLICIT  CONVOLUTION  RESULT 


ivw;  -o  - 


DIMENSION  V<200>»G<200>»  Q<200> 
S=0. 

0=0.0 

G<1>=0<1>/’<V<1>*DT> 

DO  984  I=2.N 
K=  1-1 

DO  985  M=1»K 
Q=Q+ V <K+2-M  >*G  <M> 

985  CONTINUE 
S=Q 

G<I>=<<0<I>/DT>-S>/V<1) 


Q=0.  0 

S=0.0 

984  CONTINUE 
RETURN 
END 


SUBROUTINE  EXACT1 <MM>  DT > EQVK > EQMASS > T > GG. EGVC> 


THIS  SUBROUTINE  CALCULATES  THE  REDUCED  TRANSFER  FUNCTION 
BY  DIRECT  SUBSTITUTION 


INPUT- 


EQVK--EQU I VALENT  K PARAMETER 
EQMASS— EQUIVALENT  MASS 
EQVC— EQUIVALENT  C PARAMETER 
DT— TIME  INCREMENT 
MM— NUMBER  OF  DATA  POINTS 


OUTPUT-  GG— MM  VALUES  OF  THE  TIME  DEPENDENT  TRANSFER  FUN 


DIMENSION  T<200)»GG<200) 
TRACE  ARITHMETIC 


DO  630  1=1. MM 
T< I >=FLOAT  < I >*DT 
TIME=T< I) 


BB=SQRT  < EQVK^EQMASS > 
AA=EQVC^<2. *EQMASS> 


TR I =BB*BB-AA*AA 
IFCTRI.LT.  0.  0>GO  TO  3 


GE 1 = 1 . •'SORT  < TR  I > 

GE2= 1 . /EXP  < AA*T  < I ) ) 
G3=SIN<SQRT<TRI>#T<I>> 
GG< I )=GE1*GE2*G3 


GO  TO  630 


CONTINUE 

PP=SQRT  < ABS<  TR I ) ) 
P1=EXP<PP*TCI>> 
P2=1./'<EXP<PP*T<I))) 
P3=l . /'EXP<AA*T  < I) ) 
P5=<Pl-P2>/2. 

GG<  I >=P3*P5/'PP 


CONTINUE 


RETURN 

END 


NEW  RECORD 


m 


CO  CO  CO  <0  CO  CO  o 
o o o o o o 
o o o o o o 
o o o oo  o 
o o o o o o 

r r r r r r 


cu 


o 

c 

o 


O *H  «H  OJ  (U  fU  O 
o o o o o o 
o o o o o o 
o o o o o o 
o o o o o o 

f \ r r r r 


s£  CO  CO  CU O 


o o o o o o 
o o oo  o o 
o o o o o o 


CO 

g 

o 

© 


a 


3 

a! 


h 

s 


o 

o 

♦ 


(0 


t 

w 

s 

cu 

o 

• 

S0 

t 


a 

a 


*H  *-4  O O O 
o o o o o 
I I + *f  -f 

onooo 
0 (0  o o c 
o 00  o o o 

to  eo  io  to  in 

• ••0m 

<um-*cuco 


i 

to 

<r 

£ 


C0£>  CO  ^ 
o o o o o 
I I I I I 

uuiilute 

«-<  m in  cu  co 

ID  f>w  0-t  Tt  ^ 

locrvcoocu 

C0*“  sOCO^J 

*■  co  ^ rj 


CT'  ID  CO  CO  CO  W O 

2*-  *“*“*- *• 
o o o o o o 
o o o o o o 
o o o o o o 


o 

CO 

o 

> 

COCO  (0  *“  *- 
ooooo 

in 

Ct 

i 

1 t 1 1 1 

CO 

UI 

UI 

Q 

UJ  UI  UJ  UI  Ui 

o 

b- 

*■ 

UI 

|">  ^ pyj  g 

O 

y 

sO 

»- 

0 

<0 

<r 

lOCOWNCO 

<T 

CT* 

j 

to  \0  — CO  (V 

% 

r* 

3 

CO  lO  N CO  O 

<r 

• 

u 

• • • • ■ 

a. 

'f  CO  -<  cC  CU 

<r 

8 

(0  CU  CTk  CTm»  so  O 

so  so  in  in  in  id 
o o o o o o 
o o o o o o 
o o o o o o 


•4 

CO 

*■ 

o 

o 


233 


suw  r>  v:.%  -f*  'y^  i c7  - r *■ ' riTFi-r.  tvr-.tr  n- 


■ VKimii  ujri 


V 'fc  ' 


M oooooo  oooooooooooooooooooooo 


o oooo 


> v 

m 

k£ 


SaxyyjCTi^osio^-o^  — NNoo  — ioojNNo  — — 
vpoNvDcrxoeo^^'i-ooo^^oco^Tmio  — m'Xioo^ooio 

ff.OOOO^<£i<ViOO»<N<\Ji5a)>i)N-<IO  — COtf  <UN<T>0  — tinn 
CO0MnNCO*,*-*<M'-min(M'-^'OJO''OOff'(U<0'f)CO 

N«oo<»N<jNioi^y5oajojir)>fO'i5if5^-<jS'$ooj(jsir5NnajoN 
-n-(D(con  N-Hincr.00 100 — in^^-N'-cy 
vO'CvDI^.OM-ocoosOOONffivOOIONCJNIoaiNvDaivD^’C-^vD'i 
•I  «H  (U  (U  <0 m y)  00  O (U  CO  sD  00  V ^ CO  •<  N »■•  Tf  ’-•  CO 


>i^*-*-ogo*4rcr>ajootf©c')anroMnc')co<\moo«-*wcii 
><uoo  — (yoinn^ON.-xuin’fr  — aoajoooo  — ocuoooo  — m© 

> *-  <u  <o  in  n eooNnjono«-M.NO(V4-(Di£i(iinc(uo>c 

KUOSm^tftWNOON 

)-H(u«M#ioninNo 


OM^ocy^-oo®^- 
00vDsDs0r-vD*-<CU00IO5f0'N 
00  00  00  <” 


■a-moo  — aHjNinoo^-Tj-^-oosoi^voininina) 
oo Is- no  o o»o»iooi'-nnoo(y^'Tr8soo>n 
If)  I'-  vD  *0  os  <\J  vD  7*'J0C0lf)00n00<VJO|v.|f)N0Sv0 


O»M0nj-IDNO( 

(Ms!  (u  ^ ai 


>o 

><U'f)y«'JOmOOOO)OOaiONIONO)vOWOO(<)N 

»-<«HW<vin^Mr)^CT'  — oovNoMy^-ooftjoovio 

Ti  H *H  H H OJ  fU  (\J  (\J  *H  rl 


If)  SloaNN 

SwjO 


'Dos'0'Oos(\iif)ftiio^oo'f)NvoioooN(u^-aioo^-nooo)<ugoi^<u 
v ■*■  ■*  •£  -•oofoai^ai'SooNnooi^oOMv-^-oco  — S^-oo 

^N^^r^.iiaoiOI^OOOOOIV-^-l^^-ocolO'JJOOvD'C^-O'aKMf) 

8ivioif)i^o(goooio^«fl^-<j)M^ir)vO'OU'^voo«^ir)^-  — oon 
oowwooin'tfflooAii^cu^ooajN  — — <3n°vo^-nn»-«^nn 

(UC)CU<vj<,)^l0^0)>-<lf)<=00^00<yv0l00)0)(U0)00^'jNNO0) 

^•'5-*^-'D®^ilf)l^(UO(M000lf0^-O05NI^  — MOOOCO'.XI(U<f>N 
■H  — — ainooin'ONOono^n^-  — o)^--*n^-*oo 


| ri  H H 

I I I 


^„„„w<ucuan-'fioo<r>'Oinir>cooocr>oo'0inaj®<DOinco*- 


CU  ’i 

*-ooa> 


(u^-^''Dooooaxu«i^-anoooNooooioo)nioN®N<U'£iiO'£iioo 

i'-aiaioo*"i-r'iM^-«-o'"£-<-<>0'ri'floO"H'jo-«o^-(jsaiN.ovo<\j 

— cuooai  — — ajcunwajioNoonooinai  — 

i i i i i i i 

i i 


■^(U00^-I0>f)l^'00  0)0  — ai00^-|0>0l^000)O-4(U00^l0'£N00<T' 


235 


mu 

Li  -J>  O,  -o 

1. 


PROGRAM  PRF.DCT  < I MPUT» OUTPUT » TAPE3-INPUT.  TAPE6-0UTPU1'.  PLOT  > 


2. 

3. 

4. 
3. 
6. 

7. 

8. 

9. 

10. 

11. 


12. 

13. 

14. 

13. 

16. 

17. 

18. 
19. 


20. 

21. 

22. 

23. 

24. 
23. 


C 


C THIS  PROGRAM  PREDICTS  THE  SIGNATURE  AT  AN  UNTESTED  SITE 

C DUE  TO  VEHICLE  A. IF  THE  SIGNATURE  OF  A STANDARD  VEHICLE 

C IS  KNOWN  AT  THE  UNTESTED  SITE.  AND  THE  SIGNATURE  OF  BOTH 

C THE  STANDARD  VEHICLE  AND  VEHICLE  A ARE  KNOWN  AT  A 

C STANDARD  SITE 

C 

c 

C STANDARD  VEHICLE  AND  SITE 

C C 131  SITE  1 47  SIGNATURE 

C 

C 

C C 131  SITE  3 47  BOYER  SIGNATURE 

C C-135  SITE  1 47  BOYER  SIGNATURE 

C 

C 

C INCLUDED  ONLY  FDR  PURPOSES  OF  THIS  EXAMPLE-THIS  SIGNATURE 

C IS  THE  PREDICTED  SIGNATURE 

C C 133  SITE  3 BOYER  SIGNATURE  47 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


INPUT- 

SCRO 

SIGNATURE  DUE  TO  STANDARD  VEHICLE 
AT  STANDARD  SITE 

X<  ) 

SIGNATURE  DUE  TO  STANDARD  VEHICLE 
AT  OTHER  SITE 

DFCE<  ) 

SIGNATURE  DUE  TO  OTHER  VEHICLE 

AT  STANDARD  SITE 

CDMPV<  ) 

SIGNATURE  DUE  TO  OTHER  VEHICLE 

AT  OTHER  SITE<FOR  PURPOSES  OF 

THIS  EXAMPLE  ONLY) 

ALL  INPUT 

SIGNATURES  ARE  IN  INCHES 

OUTPUT 

DO<  ) 

PREDICTED  SIGNATURE  IN  FT. 

SUBROUTINES 

REQUIRED 

<1>  EXACT 1 < ) 

<2>  NWCDNV<  > 

<3>  ERECT2<  ) 

<4)  MODPRM<  ) 

DIMENSION  SDUTK200).  SOUT2<200)>VDUT1(200) 

DIMENSION  ETINP<200),EQFNC<200) 

DIMENSION  ET <200) 

DIMENSION  VISIG<200) 

DIMENSION  XT<6» 200) . XF<6>  200) . DTRK1 <200) . DF4<200> . SFNC<200> 
DIMENSION  PP<800).XG<800).00<800) 

DIMENSION  DRV0< 1 00). DRV1 < 100>»DRV2< 1 00) 

DIMENSION  DF(200).X(100).DFCE<100).DDEF(100) 

DIMENSION  SCR(100).T<200). V(200).D<200).G(200).GG(200> 

DIMENSION  CDMPV<200).VOUT2<200) 

C*  TRACE  ARITHMETIC 

C 

C 

30  FORMA T< ^6 0X» 2 OH ") 

498  FORMAT < 29X . 34HEQU I V INPUT  IN  LBS  TDT  FUNCTION  RESPONSE  IN 
iFTss) 

398  F0RMAT<31Xi£15.3.5X.£15.3>  2X.E13.5) 

C 

READ<3> 1 00) <SCR< I ) » 1*1. 47) 

READ<3. 1 00) <X< I ) . 1*1 . 47) 

READ<3» 100) <DFCE< I ). 1*1.4?) 

READ<5. 100)<COMPV<  I).  1-1.47) 

100  FORMAT<7F10.6> 

C 

C 


DO  496  1*1.47 
SOUT1 < I )*SCR< I >•<'12. 
V0UT1<1)*DFCE<I)/12. 
S0UT2<I)«X<I)/12. 
VOUT2<I )-COMPV< I)/12. 
496  CONTINUE 
C 


JL 


237 


UAJJ 


li .O'. 


i»; 

ii 


-j  w-jjtts.  w.'J 


rw  rrjrraw  rar  rv*  pl»  rv*  rr»  r » rrr  n v»n 


TTW  k 


JP=47 

PMAX  =20000. 

DT=0. 04 
MM=47 
LL=47 
NT =47 

DO  74  1=1 i NT 
VISIG<I>=S0UT1<I> 
74  CONTINUE 


WRITE<6,497> 

497  FORMAT <.//£Xr  23HS I TE  ONE  STD  VEHICLE/'/') 
CALL  MODPRM  <PMAX, JP , DT  > V I S I G , EQMASS , DF  > GG ) 


COMPUTE  EQUIV.  INPUT  FOR  VEHICLE  2 SITE  1 
DO  721  1=1, MM 
T <1 )=FLQAT  < I )*DT 
□ <I)=VOUTia> 

V<I)=GG<I) 

721  CONTINUE 


N=47 


CALL  NWCONV<NiViG»O.DT) 

DO  722  1=1, MM 
ETINP< I )=G< I )#EQMASS 
722  CONTINUE 

URITE<6,50) 

WRITE<6, 495) 

495  F0RMAT</'/'2X,21HSITE  ONE  VEHICLE  2/-^) 
14RITE<6, 498) 

14RITE<6,398)<ETINP<I)»GG<I)»Q<I>» 1=1, MM) 


COMPUTE  EQUIVALENCY  FUNCTION 
DO  723  1=1, MM 
EQFNC<I  )=ETINP<  I )/”DF<  1) 

723  CONTINUE 

COMPUTE  TDT  FUNCTION  , EQUIV.  INPUT  FOR  NEXTSITE  STD.  VEHICLE 
I4RITE<6,397) 

397  FORMAT </'/'2X» 23HNEXT  SITE  STD  VEHICLE/'/') 


DO  724  1=1, MM 
VISIG<I)=SQUT2<I) 

724  CONTINUE 

CALL  MODPRM<PMAX,  JP,  DT, VIS IG, EQMASS , DF , GG) 


COMPUTE  EQUIV.  INPUT  NEXT  SITE  VEHICLE  2 

DO  725  1=1, MM 

ET < I ) =DF <1 ) *EQFNC <1 > 

725  CONTINUE 

COMPUTE  OUTPUT  FOR  NEXT  SITE  VEHICLE  2 

DO  206  1=1, MM 

PP  < I ) =ET  < I ) /'EQMASS 

XG<I)=GG<I) 

206  CONTINUE 
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69.  CALL  ERECTECMMjPP.XG.DTiOO) 

C 

C 

70.  WRITE<6*  50> 

71.  WRITEt6»494> 

7S.  494  FORMAT <//'£X»21HNEXT  SITE  VEHICLE  2/v> 

73.  WRITE<6>498> 

74.  WRITE  <6»  398  XET  (I)>XG<I)>00(I>» I«1»MM> 

C 

75.  WRITE<6>50> 

76.  WRITE<6»2000> 

77.  2000  FORMAT  < /V6 OX » »PRED ICTED*>  1 0X«  *MEASURED*/V> 

78.  WRITE<6.1000><00<:i>»VOUT2< I >.1=1.47) 

79.  1000  FORMAT <52X* 2F15.-5) 

C 

80.  WRTTE<6.50> 

C 

C 

C 

C 

C 

C 

81 . STOP 

82.  END 
C 

C 

C 

1.  SUBROUTINE  NWCONVCNj V»G»0»DT> 

2.  DIMENSION  V<200> > G<200> » 0<200> 

3.  S=0. 

4.  <2=0. 0 
C 

C 

5.  G<1)=0<  1)/’<V<  1)*DT) 

6.  DO  984  I=2»N 

7.  K=I-1 

8.  DO  985  M=1»K 

9.  Q=Q+V<K+2-M2*G<M) 

10.  985  CONTINUE 
C 

11.  S=<2 

12.  G<I>=<<0<I  VDT)-SW<1> 

• C 

C 

13.  <2=0.0 

14.  S=0. 0 

15.  984  CONTINUE 

16.  RETURN 

17.  END 
C 

C 

C 

1.  SUBROUTINE  ERECT2 < MM . PP , XG > DT » 00 > 

2.  DIMENSION  PP<800XXG<800>, 00(800) 

C 

C 

3.  PS2=0, 0 

4.  DO  331  I=2»MM 

5.  K=I-1 

6.  DO  332  M=1 >K 

7.  PS2=PP(K+1-M)*XG<M>+PS2 

8.  332  CONTINUE 

9.  OOCI-l >=PS2*DT 

10.  PS2=0. 0 

11.  331  CONTINUE 

12.  RETURN 

13.  END 
C 

C 

v 

1.  SUBROUTINE!  EXACT  ( MM » DT » EQVK » EQMASS » T » GG > EQVC ) 

2. ^*  DIMENSION  T(200)rGG(200> 
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C COMPUTING  EXACT  TRANSFER  FUNCTION 

DO  630  1=1, MM 
T< I >=FLOAT  < I >*DT 
TIME=T  < I> 

C 

C 

BB=SQRT  < EQVK/'EQM  ASS  > 

AA=EQVC^<£. *EQMASS> 

C 

TR I =BB*BB-AA#A A 
IF <TRI.LT. 0. 0>GO  TO  3 
C 

GE1=1./SQRT<TRD 

GE£=1./'EXP<AA*T<I>> 

G3=S IN  < SORT  < TR I > *T  < I > > 

GG< I >=GE1*GE£*G3 
C 
C 

GO  TO  630 
C 

3 CONTINUE 

PP=SQRT  < ABS  < TR I > ) 

P1=£XP<PP#T  < I > > 

P£=l . /’<EXP<PP*T  < I > > > 

P3=l . ^EXP<AA*T < I>  > 

P5=<P1-P£V£. 

GG< 1 )=P3*P5/PP 
C 

630  CONTINUE 
C 
C 

RETURN 

END 


SUBROUTINE  MODPRM<PMAX, JP, DT, VISIG, EQMASS, DF» GG> 

DIMENSION  QC<100> 

DIMENSION  DF <£00) 

DIMENSION  T <£00) 

DIMENSION  0<£00>,GG<£00> 

DIMENSION  V<£00) » G<£00) 

DIMENSION  IVEL< 1 1,5) 

DIMENSION  DRV1 <£00)»DRV£<£00>,  CPMXX<50>» VISIG<£00>»DRV0<£00> 
C 
C 

DO  74  1=1, JP 
DRV0< I )=VISIG< I ) 

TRM1=<VISIG<I>+VISIG< 1+1 ))^£. 

IF< I . GT. 1 )GO  TO  744 
TRM£=VISIG<IV£. 

GO  TO  474 
744  CONTINUE 

TRM£=<VISIG<I)+VISIG<I-1)V£. 

474  CONTINUE 

DR V 1 <1 > = < TRM 1 ~TRM£  > -"DT 

74  CONTINUE 

DO  75  1=1, JP 
TRN1=<DRV1 < 1+1 )+DRVl< I ) 

IF< I.GT. 1 >GO  TO  755 
TRN£=DRV  1 < I ) /'£ . 

GO  TO  575 
755  CONTINUE 

TRN£=<DRV1<I)+DRV1<I-1)V£. 

575  CONTINUE 

DRV£<  I > = <TRN1-TRN£>/'DT 

75  CONTINUE 
C 

c 

BOTM=0. 00000000001 

I;************* 

DO  7 £6  1=1, JP 

HUMP 1 =DR V 0 < I + 1 > -DRV  0< I > 

IFCHUMPl.GE. 0. 0>GO  TO  7£6 
JPP=I+6  
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33.  GO  TO *72? 

36.  726  CONTINUE 

37.  727  CONTINUE 

C************ 

c 

38.  IMAXV-0 

38.  DV1MX=0. 0 

40.  DO  313  I-l.JPP 

41.  IF<DRV1<I).EQ.  0.  0>DRV1<I>-<DRVKI+1>+DRV1<I-1>V2. 

CAUTION  - IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  ****************************> 

42.  IF<DRV1<I).EQ.  0.  0)DRVl<I)-<DRVl<I+2)+DRVl<I-2>>.'2. 

CAUTION  - IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  ***************************** 

43.  IF<ABS<DRVl<I)).LT.B0TH>DRVKI>»<DRVl<I+l>+DRVKI-l»/2. 

44.  IF<ABS<DRV2<  I >>.LT.BQTM>DRV2<  I >"<DRV2<  I +1>+DRV2<  I -DJ/'E. 

43.  IF<DRV1<I).LT.DV1MX)G0  TO  313 

46.  DV1MX=DRV1 < I > 

47.  INAXV*I 

48.  313  CONTINUE 
C 

C 

49.  DMAX=DRV0<1> 

50.  DO  33  1=2. JPP 

51.  IF<DRVO<I).LT.DMAX)GO  TO  53 

52.  DMAX=DRVO<I> 

53.  IMAXD=I 

54.  33  CONTINUE 
C 

C 

55.  IMAXV1=INAXV-1 

56.  DO  57  I=2» IMAXV1 

57.  IF<DRV2< I ) . LE. 0. 0>DRV2< I )=<DRV2< 1-1 >+DRV2CI+l > 

58.  IF<DRV2<I).LE. 0. 0>DRV2<I>=<DRV2<I-2>+DRV2<I+2>>^2. 

59.  57  CONTINUE 

60.  I MAXVA= I MAXV+ 1 

61.  DO  109  I=IMAXVA» JPP 

62.  IF  <DRV2< I ) .GE. 0. 0)DRV2< I )=<DRV2< 1-1 >+DRV2< 1 + 1 > >✓2. 

63 . I F<  DRV2  < I > . GE . 0 . 0 > DRV2 < I > « < DR V2  < I -2) +DRV2  <I+2>)^2. 

64.  109  CONTINUE 

65.  DO  107  1=1. JPP 

66.  IF^DRVl < 1+1 ) . EQ. DRV1 < I ))DRV1 < 1+1 >=<DRV1 < I+2>+DRVl < I > >>'2. 

CAUTION  - IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  **************************** 

67.  IF<DRV2<I+l).EQ.DRV2<I))DRV2<I+l)=<DRV2<I+2)+DRV2<I>>-'2. 

CAUTION  - IS  FLOATING  POINT  EQUALITY  TO  BE  EXPECTED  ****************************. 

68.  107  CONTINUE 
C**M***K 

69.  DO  789  IA=1 > JPP 

70.  DEROT  S=DRV  1 < I A+ 1 > >"DR  V 1 < I A ) 

71.  IF<DEROTS.GT. 0. 0)GO  TO  789 

72.  DZERD=DRV 1 < I A ) -DRV 1 < I A+ 1 ) 

73.  DDTZ=DRV1<IA)>-D2ER0 

74.  NAFNL=IA 

75.  GO  TO  788 

76.  789  CONTINUE 

77.  788  CONTINUE 

C««**MM**** 

78.  DO  889  IA=t>  JPP 

79.  ZER0TS=DRV2< IA+1)/'DRV£< IA> 

80.  IF<ZEROTS.GT. 0. 0>GO  TO  889 

81.  TZER0=DRV2< I A)-DRV2< IA+ 1 > 

82.  FDTZ=DRV2<IAVTZER0 

83.  IAFNL=IA 

84.  DR0RNG=DRV0<IA+1)-DRV0<IA> 

85.  ZVALU  0=DRV  0 < I A > +FDTZ*DR  0RNG 

86.  DRlRNG=DRVl<IA+l>-DRVl<IA-> 

87.  ZVALU 1=DRV1 < IA)+FDTZ*DR1RNG 

88.  GO  TO  888 

89.  889  CONTINUE 

90.  888  CONTINUE 

C 

C**««****» 

c 

c 

91.  WRITE<6i 99) 

92.  99  FORMAT < ^30X > *FOR  SIGNATURE*/') 

93.  WRITE<6»  76) 


i 

■ i&m&rk  mm 


lT»%srru  ■ »TT»H  « u ■ u ■ 1 HI!  UM  HJLWJlUJt  M U.miMUBaanMi  u>  I WIIU.W1BIB  I u>»  _■  ijj.ji  ) 


94.  76  FORMAT</vlX*“  INCREMENT  DEFLECTION  VELOCITY 

1ERATIQN*<0 

95.  URITE<6»77>< I>DRV0< I>>  DRV1 <I>»  DRV2< I>» !■!» JPP> 

96.  77  F0RMATC1X. I14i3E15.5> 

97.  78  F0RMAT</'15X»2I15/'> 

Cmxxxxxxxxxxixxxxi 

CKXWXIXUKXMXI 

98.  IDIFF=IMAXD-IMAXV 

99.  IMAXF=0 

100.  AA1«=0. 05 

101.  ICC»0 

CKMMMXnXNNKXl 

102.  DO  104  ILP~1»40 

C*****MMM 

103.  EQKA=AA1*PMAX/'DMAX 

104.  AA3--0.91 
0*** 

105.  DO  103  JJ*1»20 

106.  IT*IMAXV 

107.  EQMA=AA3“PMAX/'DRV2<  IMAXD) 

108.  CPMAX=0. 0 

109.  AA2=AA1+AA3 

110.  IF<AA£.LT. 0.0>GO  TO  1021 

CKXMM 

111.  IMAXC=IMAXD 

112.  IT=IT+1 

113.  DO  102  I I=IT» IMAXC 

114.  EQCA=EQMA*DRV2 < I 1 VDRV1  < 1 1 > 

115.  EQCA=-EQCA 
C 

C 

116.  EQVK=EQKA 

117.  EQMASS=EQMA 

118.  EQVC=EG(CA 

1 19.  MM= JPP 

120.  CALL  EXACT 1 < MM » DT » EQVK » EQMASS > T i GG i EQVC > 

C 

C 

C 

121.  DO  73  I=1»JPP 

122.  0<I)=DRVO<I> 

123.  V<I)=GG<I) 

124.  73  CONTINUE 
C 

125.  N=JPP 

126.  CALL  NWC0NV<N.V.G»0.DT) 

C 

127.  DO  508  MP=1 » JPP 

128.  DF<MP)“G<MP> “EQMASS 

129.  508  CONTINUE 

130.  DO  509  MA= 1 * JPP 

131.  DF  <MA)=<DF  <MA>+DF<MA-*1 ) )^2. 

132.  509  CONTINUE 
C 

C 

133.  CPMAX=0. 

134.  DO  501  M=1 » JPP 

135.  DF  <M>=<DF<M>+DF  <M+1>>^2. 

136.  IF<DF<M>.GT.CPMAX)GO  TO  20 

137.  GO  TO  502 

138.  20  CPMAX=DF<M> 

139.  ICF=M 

140.  502  CONTINUE 

141.  501  CONTINUE 
C 

142.  ERFCE=CPMAX-PMAX 

143.  AEERFG=ABS  < ERFCE  > 

144.  IF<ABERFG.GT.500. 0>GO  TO  102 
C 

145.  ICC=ICC+1 
C 

146.  DFRNG=DF ( IAFNL+1 )-DF  < IAFNL) 

147.  2VALUF=DF< IAFNL )hFDT2“DFRNG 

1 48 . CMAXV= < 2VALUF -EQKA*2VALU 0 > /2VALU1 

149.  CD I FFR= (CMAXV-EQCA ) /"CMAXV 


u*  m!  twi  lp i in  j.m n w mi  u»  i m i iwi  iw  i n 
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150. 

ABCD=CDIFFR 

151. 

QC< ICC>=CDIFFR 

152. 

IF < ICC. EQ. 1 >GO  TQ  102  . 

153. 

CRATIO=QC<  ICC>/'QCaCC-l  > 

154. 

C 

r- 

IFCCRATIO.GT.O. 0>GO  TQ  102 

155. 

L > 

c 

FCHK=EQKA*DRV0< I I > 

156. 

DIFCHK=FCHK-DF  < 1 1 > 

15?. 

DFCHKR=D I FCHK/FCHK 

158. 

FPFRAC=DF<II>/'PMAX 

159. 

ICDPT=I I 

160. 

I BBB= I MAXD- I CDPT 

161. 

BB1=FL0AT  < IBBB^IDIFF) 

168. 

TRMM=EQMA*DRV2< I I > 

163. 

TRMC=EQCA*DRV1 < I I > 

164. 

TRMK=EQKA*DRV0< I I > 

165. 

CPMXX< 1 1 > =TRMM+TRMC+TRMK 

166. 

WRITE<6»847> 

16?. 

84? 

FORMAT < /v  1 0X>  *G00D  RUN*/V> 

168. 

GD  TQ  846 

169. 

102 

CONTINUE 

170. 

1021 

CONTINUE 

171. 

AA3=AA3+0. 045 

QK*K»K«*K*K *** 

172. 

103 

CONTINUE 

173. 

AAl=AAl+0. 025 

174. 

104 

CONTINUE 

175. 

1000 

FDRMAT</'1X»*000*> 

176. 

844 

FDRMAT<1Xj7F10.2»3E15.5> 

177. 

845 

F0RMAT<6F 15. 4j 10X»2I5.F10.3> 

178. 

URITE<6> 1002) 

179. 

1002 

FORMAT < ^ 1 OX j *OUTER  LOOP  COMPLETED*^) 

180. 

C 

GO  TO  94 

181. 

c 

846 

r- 

CONTINUE 

182. 

URITE<6j 1000) 

183. 

AFRNG=DF  <NAFNL+1 >-DF<NAFNL> 

184. 

DVALUF =DF  < NAFNL  ) +DDTZ*AFRNG 

C«*KKK 

185. 

ZVALUR=ZVALUF/PMAX 

C**KKKK 

186. 

DVALZD=DVALUF'PMAX 

£K*KKK 

187. 

FDRT IO=ABERFG'PMAX 

Qmmmmmmmm 

188. 

DENOM=FLOAT < 1MAXD> 

189. 

DENOMl=l. 

190. 

TCRSPT=FLOAT  < 1 1 VDENOM 

191. 

TZROA=FLOAT < I AFNL  > 'DENOM 

CWKKKKKK  * 

192. 

1199 

FORMAT</'10X.  11F11.3^> 

193. 

1001 

FORMAT < 1 0X» I5»  7F15. 5> 

194. 

C ALL  EXACT 1 < JP > DT  > EQKA  t EQMA » T > GG i EQC A > 

195. 

DO  733  1=1 »JP 

196. 

0<I>=DRVO<I> 

197. 

V<I)=GG<I) 

198. 

733 

CONTINUE 

199. 

N=JP 

200. 

CALL  NWCONVCN*  V>  Gi 0»  DT> 

201. 

DO  518  MP=1 > JP 

202. 

DF<MP>=G<MP>*EQMA 

203. 

518 

CONTINUE 

204. 

C 

205. 

94 

CONTINUE 

206. 

RETURN 

207. 

END 
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SITE  ONE  STD  VEHICLE 


FOR  SIGNATURE 


INCREMENT 

DEFLECTION 

VELOCITY 

ACCELERATION 

1 

1 . 27667E-04 

1.92708E-03 

9.93490E-03 

2 

1 . 54167E-04 

7.  94792E-04 

8. 31380E-03 

3 

1 . 91250E-04 

1. 03125E-03 

6. 69271E-03 

4 

2 . 36667E- 04 

1. 33021E-03 

1. 03906E-02 

5 

2. 97667E-04 

1. 8625 0E- 03 

1. 57552E-02 

6 

3. 85667E-04 

2. 59063E-03 

2. 1 0677E-02 

7 

5.  04917E-04 

3. 54792E-03 

2. 91276E-02 

8 

6. 69500E-04 

4. 92083E-03 

3. 62500E-02 

9 

8. 98583E-04 

6. 44792E-03 

4. 09245E-02 

10 

1. 18533E-03 

8. 19479E-03 

4. 38151E-02 

11 

1.55417E-03 

9. 95312E-03 

3. 51432E-02 

12 

1.98158E-03 

1. 10062E-02 

9. 81771E-03 

13 

2. 43467E-03 

1. 07385E-02 

-3. 29948E-02 

14 

2.84067E-03 

8.  36667E-03 

-8.32943E-02 

15 

3. 1 0400E-03 

4. 07500E-03 

-1. 15703E-01 

16 

3.  16667E-03 

-8. 89583E-04 

-1. 14401E-01 

17 

3. 03283E-03 

-5. 07708E-03 

-8. 50130E-02 

18 

2.  76 05 0E- 03 

-7. 69062E-03 

-4. 40234E-02 

19 

2. 41758E-03 

-8. 59896E-03 

-4.75260E-03 

20 

2. 07258E-03 

-8. 07083E-03 

-9. 53125E-03 

21 

1.77192E-03 

-6. 89896E-03 

7. 39583E-03 

22 

1 . 52067E-03 

-5.80104E-03 

2. 50000E-03 

GOOD  RUN 


SITE  ONE  VEHICLE  2 


> INPUT  IN  LBS 

TDT  FUNCTION 

RESPONSE  IN  FT 

7. 26622E+03 

2. 30853E-02 

2. 02000E-04 

9. 67190E+02 

2. 76124E-02 

2. 68500E-04 

3. 38946E+03 

2. 56125E-02 

3. 50500E-04 

6. 28019E+03 

2. 17664E-02 

5. 07583E-04 

1. 20998E+04 

1 . 78084E-02 

8.30917E-04 

2. 26201E+04 

1 . 43096E-02 

1.45967E-03 

3. 73522E+04 

1. 13959E-02 

2. 51742E-03 

5. 36091E+04 

9. 03463E-03 

4.03258E-03 

6. 66077E+04 

7. 1461 IE-03 

5. 86650E-03 

5. 49146E+04 

5.64567E-03 

7. 24833E-03 

-4.4064 0E+03 

4. 45756E-03 

6.66300E-03 

-1. 88870E+04 

3. 51838E-03 

5. 11650E-03 

1 . 76576E+04 

2. 77663E-03 

4. 47067E-03 

4. 69750E+04 

2. 19108E-03 

5. 00992E- 03 

8.37103E+04 

1.72893E-03 

6 . 88358E- 03 

6. 44724E+04 

1 . 36424E-03 

8. 41667E-03 

-1 . 57651E+04 

1 . 07645E-03 

7.41800E-03 

-2. 57654E+04 

8. 49373E-04 

5.45308E-03 

-5. 79904E+03 

6. 70194E-04 

3. 97867E-03 

-1.63126E+03 

5. 28813E-04 

2.96208E-03 

1 . 61962E+03 

4. 1725/E-04 

2. 31 008E-03 

4. 38986E+03 

3.29234E-04 

1.93375E-03 

4. 45170E+03 

2. 5978 0E- 04 

1.69475E-03 
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5.24889E+03 
5. 77634E+03 
4.82121E+03 
5. 39067E+03 
6. 05333E+03 
4. 47763E+03 
4.63410E+03 
5. 15415E+03 
5. 09242E+03 
3. 93199E+03 
5. 04313E+03 
5. 26145E+03 
4. 01299E+03 
4. 74872E+03 
6. 32678E+03 
2.59107E+03 
5. 12136E+03 
5. 51294E+03 
4.25295E+03 
4.31075E+03 
5. 95350E+03 
4. 67396E+03 
3. 94255E+03 
4. 68328E+03 


2. 04977E-04 
1.61736E-04 
1.27617E-04 
1. 00695E-04 
7.94527E-05 
6. 26916E-05 
4.94664E- 05 
3.90311E-05 
3. 07972E-05 
2.43003E-05 
1.91740E-05 
1 . 51291E-05 
K19375E-05 
9.41922E-06 
7.43217E-06 
5.86430E-06 
4.62719E-06 
3. 651 05E-06 
2.88083E-06 
2. 2731 0E- 06 
1.79358E-06 
1 . 4 1 52 1 E- 06 
1. 11666E-06 
8.81 093E-07 


1.55192E-03 
1 . 4723 0E- 03 
1 . 39683E-03 
1.34767E-03 
1 . 33158E-03 
1.28433E-03 
1.23733E-03 
1.21075E-03 
1. 19233E-03 
1 . 14658E-03 
1 . 12867E-03 
1. 12800E-03 
1 . 09825E-03 
1 . 08333E-03 
1. 11892E-03 
1 . 06242E-03 
1 . 05375E-03 
1 . 07242E-03 
1. 06250E-03 
1 . 04625E-03 
1 . 07567E-03 
1. 08050E-03 
1. 0565  0E- 03 
1.04683E-03 


S! 


NEXT  SITE  STD  VEHICLE 

FOR  SIGNATURE 


INCREMENT 

DEFLECTION 

VELOCITY 

ACCELERATION 

1 

9.95833E-05 

1.52187E-03 

7. 77344E-03 

2 

1.21750E-04 

6.21875E-04 

6.71224E-03 

3 

1.49333E-04 

8. 09375E- 04 

5.65104E-03 

4 

1 . 86500E-04 

1.07396E-03 

6. 19792E-03 

5 

2.35C50E-04 

1.30521E-03 

5. 92448E-03 

6 

2.90917E-04 

1.54792E-03 

7. 1224  0E- 03 

7 

3.59083E-04 

1.87500E-03 

8.54167E-03 

8 

4.40917E-04 

2.23125E-03 

9. 60937E- 03 

9 

•5.37583E-04 

2.64375E-03 

1.25000E-02 

10 

6. 52417E-04 

3.23125E-03 

1 • 47917E- 02 

11 

7.96083E-04 

3.82708E-03 

1.32161E-02 

12 

9.58583E-04 

4.28854E-03 

7.79948E-03 

13 

1. 13917E-03 

4. 451 04E-03 

-6. 17187E-03 

14 

1.31467E-03 

3. 79479E-03 

-2. 66797E-02 

15 

1.44275E-03 

2.31667E-03 

-4 . 29948E-02 

16 

1.50000E-03 

3. 55208E-04 

-4.88021E-02 

17 

1.47117E-03 

-1. 5875 0E- 03 

-3. 8 099 0E- 02 

18 

1.37300E-03 

-2.69271E-03 

-1.61849E-02 

19 

1.25575E-03 

-2. 88229E-03 

-1. 82292E-03 

20 

1. 14242E-03 

-2. 83854E-03 

-5. 88542E-03 

21 

1 . 02867E-03 

-2. 75000E-03 

-7.35677E-04 

22 

9. 22417E-04 

-2. 56146E-03 

-7. 16146E-04 

GOOD  RUh 

\S 


NEXT  SITE  VEHICLE  2 


EQUIV  INFUT  IN  LBS 

TDT  FUNCTION 

RESPONSE  IN  FT 

2. 39378E+04 

. 2. 79870E-02 

1 . 57565E-04 

-5.99763E+04 

4. 02600E-02 

-1.68119E-04 

6.45929E+03 

4. 45935E-02 

-2 . 74326E- 04 

1 . 36858E+04 

4.49856E-02 

-2. 24518E-04 

2. 55768E+04 

4. 34868E-02 

-2. 40447E-05 

3.73268E+04 

4. 1 1442E-02 

3. 17975E-04 

5. 95647E+04 

3. 84877E-02 

8. 60932E-04 

7. 51 046E+04 

3. 57813E-02 

1 . 58148E-03 

8. 03293E+04 

3.31525E-02 

2.39395E-03 

6. 33737E+04 

3. 06592E-02 

3. 1 0814E-03 

-4. 96779E+03 

£ . 83238E-02 

3. 29894E-03 

-1 . 96412E+04 

2. 6151 0E-02 

3. 13688E-03 

1.84098E+04 

2. 41370E-02 

3.  0634 0E- 03 

4. 64276E+04 

£ . 2274  IE-02 

3.21925E-03 

6. 79205E+04 

2. 05529E-02 

3. 61990E-03 

3. 74124E+04 

1 . 89636E-02 

3.921 25E- 03 

-3.1111 0E+01 

1 . 74966E-02 

3.91757E-03 

2. 12723E+04 

1 . 61429E-02 

3. 90907E-03 

5. 20213E+03 

1 . 48937E-02 

3. 79282E-03 

1. 06867E+03 

1 . 3741  IE-02 

3. 60242E-03 

-1.67619E+0? 

1.26777E-02 

3. 36579E-03 

-3. 75658E+03 

1.  16966E-02 

3. 1 0234E-03 

-4. 38155E+03 

1. 07913E-02 

2. 83187E- 03 

-6. 42824E+03 

9. 95617E-03 

2. 55472E-03 

-9. 40418E+03 

9.  18563E-03 

2. 26518E-03 

-6. 84531E+03 

8 . 47472E-03 

1 . 99745E-03 

-6. 0481 1E+03 

7. 81883E-03 

1 . 75536E-03 

-7.78484E+03 

7.21371E-03 

1 . 52313E-03 

-4. 78390E+03 

6. 6554 IE- 03 

1 . 32403E-03 

-2.34767E+03 

6. 14033E-03 

1. 16421E-03 

-3. 53550E+03 

5.6651  IE-03 

1. 02124E-03 

-8. 24185E+02 

5. 22666E-03 

9. 09508E-04 

3. 66539E+03 

4. 82215E-03 

8.46376E-04 

5. 42162E+03 

4. 44895E-03 

8.20303E-04 

-7.56887E+03 

4.  1 0463E-03 

7. 27339E-04 

7.55546E+03 

3. 78696E-03 

7.05572E-04 

1.67329E+03 

3. 49387E-03 

6.  7979 IE- 04 

4. 12036E+03 

3. 22347E-03 

6. 69173E-04 

2. 13777E+03 

2. 97399E-03 

6.53120E-04 

1.94469E+03 

2. 74383F-03 

6. 3381 0E- 04 

3. 01993E+03 

2. 53147E-03 

6. 20752E-04 

6. 73974E+03 

2. 33555E-03 

6.35642E-04 

2. 83888E+03 

2. 15480E-03 

6.37603E-04 

4. 24396E+03 

1.98303E-03 

6.42583E-04 

3. 23174E+03 

1 . 83417E-03 

6.  42152E-04 

6. 86666E+03 

1 . 69222E-03 

6.63086E-04 

2. 46450E+O3 

1.56125E-03 

0 

PREDICTED 

MEASURED 

. 00016 

. 00018 

-. 00017 

.00028 

-. 00027 

. 00044 

-. 00022 

. 00068 

-. 00002 

. 00107 

. 00032 

. 00165 

. 00086 

. 00252 

. 00158 

. 00357 

. 00239 

. 00404 

.00311 

. 00367 

. 00330 

. 00318 
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. 00314 

. 00289 

. 00306 

. 00284 

. 00322 

.00314 

.00362 

.00382 

.00392 

.00433 

.00392 

.00399 

.00391 

.00320 

.00379 

.00244 

.00360 

.00178 

. 00337 

. 00128 

. 00310 

. 00091 

.00283 

. 00067 

.00255 

. 00051 

. 00227 

.00041 

.00200 

. 00035 

. 00176 

. 00032 

. 00152 

.00029 

.00132 

. 00027 

.00116 

. 00025 

.00102 

.00014 

. 00091 

.00012 

.00085 

.00012 

.00082 

. 00012 

.00073 

.00012 

.00071 

.00012 

. 00068 

. 00G1 1 

.00067 

. 00012 

. 00065 

. 00011 

.00063 

. 00011 

. 00062 

. 00011 

. 00064 

.00011 

.00064 

. 00011 

. 00064 

. 00011 

. 00064* 

.00011 

. 00066 

.00011 

0 

.00011 

APPENDIX  D 


DATA  DIGITIZED 
FROM  LIGHT 
SENSITIVE  PAPER 

The  data  listing  has  the  following  format 

Card  1 - A comment  card  giving  the  base  name,  site  number,  grid 
location  of  measurement  point,  Load  Cart  data  of  test, 
distance  from  closest  measurement  gage  to  edge  of  tire 
track  in  inches,  time  increment  at  which  data  was  digitized 
and  estimated  number  of  cards. 

Card  2 - Distance  (ft)  of  measurement  point  from  edge  of  tire  print 
where  Gage  1 is  the  closest. 

Card  3 - Number  of  cards  actually  read,  time  increment  used,  and 
weight  of  load  cart  (25  kips). 

Card  4 to  (n-l)  - Data  digitized  and  listed  (from  left  to  right)  as. 
Gage  1,  Gage  2,  Gage  3,  Gage  4,  and  Gage  5* 

Card  n (last  card)  - Calibration  data  for  Gages  1 to  5 also  listed  from 
from  left  (Gage  l)  to  right  (Gage  5)* 
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J :.r::iN 

SITE  1 

LOG  8- 

-.G-,2 

.9375 

1.6875 

48 

m 

200  25000.0 

0027 

0045 

0071 

0102 

0134 

0030 

0048 

0074 

0106 

0134 

0034 

0051 

0078 

0110 

0136 

0040 

0058 

0085 

0117 

0137 

0043 

0065 

0094 

0126 

0141 

0061 

0078 

0107 

0138 

0145 

0077 

0092 

0122 

0154 

0151 

0103 

0115 

0146 
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